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Sommaire 


On rend compte de mesures de transparence acoustique faites sur des hublots fermés par des 
calottes sphériques de différentes courbures. Par rapport au disque plan, il y a une nette amélio- 
ration de l’isolement, mais aux graves seulement. Ce résultat est expliqué par l’étude de la répar- 
tition des fréquences propres d’une calotte vibrant en extension. 


Zusammenfassung 


Es wird iiber die Messung der Schallisolation von Kugelschalen verschiedener Kriimmung 
berichtet, die bei tiefen Frequenzen eine Verbesserung der Schallisolation gegeniiber ebenen 
Platten ergeben. Das Ergebnis wird durch eine Betrachtung der Verteilung der Eigenschwingun- 
gen der Kugelschale (bei Dehnschwingungen) erklart. 


Summary 


The acoustical insulation of spherical shells has been found to be appreciably greater than that 
of plane discs in the lower frequency range. This result is accounted for in a theoretical discussion 
of the natural frequencies of the shell, under the assumption of extensional vibrations. 


Introduction 


Les qualités bien connues des calottes sphéri- 
ques élastiques en ce qui concerne leur résistance 
a la déformation statique avaient conduit en 1939 
M. L. LumiéRreE a proposer [1] la fabrication de 
miroirs sphériques en plexiglas: ceux-ci, repo- 
sant par leur bord sur un cercle métallique par- 
faitement rigide, étaient pratiquement indéfor- 
mables. Des réalisations furent faites par l’un de 
nous en collaboration avec M. LUMIERE, et per- 
mirent de montrer qu’on pouvait, 4 condition de 
ne pas se montrer trop exigeant, remplacer le 
verre par de la matiére plastique. 

L’un de nous ayant pensé que ce procédé pou- 
vait étre intéressant pour obtenir des parois 
rigides quoique trés légéres, ne vibrant pas sous 
Veffet des ondes sonores, nous demanddmes 4 
M. LumiéreE de nous préter un de ses miroirs 
(celui désigné plus bas par le sigle L); les pre- 
miers résultats furent encourageants. Des essais 
faits par la suite sur des calottes métalliques 
montrérent l’influence capitale de la fixation du 
contour: la partie sphérique avait, en effet, été 
emboutie dans un panneau plan plus étendu, que 
Yon fixait le long de ses bords et qui vibrait lui- 
méme, d’ot pour la calotte des conditions au 
contour difficiles 4 préciser, et d’ailleurs peu 
favorables. 

On a repris récemment ces mesures en serrant 
soigneusement la partie plane entre des ma- 


choires circulaires rigides: les résultats en ont 
paru suffisamment intéressants pour mériter un 
bref compte-rendu. 


1. Description des expériences 


Les mesures ont porté sur quatre calottes, dont 
trois en aluminium et une en plexiglas, ayant 
toutes un contour circulaire de 27 cm de rayon, 
mais des courbures différentes. Les premiéres ont 
été réalisées par martelage sur un moule creux en 
ciment, et présentent des irrégularités de cour- 
bure inférieures 4 4% ; la derniére a été moulée a 
chaud et est d’une belle régularité, mais porte au 
pole un trou de 1 mm environ de diamétre. Le 
Tableau suivant indique leurs caractéristiques 
principales. 


Tableau I 
Epai R de|R du] Demi- 
Calatie Matiére hee Soebune canter darren 
7 en mm |r en cm | aen cm |@, en degrés 
A Aluminium 1 35 27 51 
B Aluminium ] 42 27 40 
Cc Aluminium 1 84 27 19 
L Plexiglas 2 79 27 20 


Serrés, comme il a été dit, sur leur contour 
entre des plaques suffisamment rigides pour qu’on 
puisse considérer comme immobiles les points du 
petit cercle de rayon a, et comme négligeable la 
transparence de la partie plane, les hublots 
étaient fixés dans la baie de mesures qui sépare 
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la grande salle sans écho du laboratoire d’une 
petite arriére-salle sourde; une source d’ondes 
planes de fréquence continuement variable et de 
niveau constant irradiait la baie sous incidence 
normale, et le son parvenu dans la chambre 
sourde, capté par un microphone, commandait 
un enregistreur de niveau Briiel et Kjzr a sensi- 
bilité logarithmique, dont le déroulement était 
asservi a la fréquence de la source. 

Une courbe de réponse analogue était ensuite 
obtenue, la calotte sphérique étant remplacée 
par un panneau plan de mémes caractéristiques 
élastiques, semblablement serré : la comparaison, 
abscisse pour abscisse, des deux courbes a permis 
de tracer les caractéristiques des Fig. la—d, qui 
mettent en lumiére l’influence de la sphéricité du 
hublot sur son isolement acoustique. 


150 300 


1000 2000 3000 4000 Hz 


7000 2000 3000 4000 Hz 


-20 


Fig. 1. (a) Isolement acoustique du hublot sphérique A 
rapporté a celui du hublot plan correspondant. 

(b) Isolement acoustique du hublot sphérique B 
rapporté a celui du hublot plan correspondant. 

(c) Isolement acoustique du hublot sphérique C 
rapporté a celui du hublot plan correspondant. 

(d) Isolement acoustique du hublot sphérique L 
rapporté a celui du hublot plan correspondant. 


On constate en gros une amélioration d’isole- 
ment aux graves, de l’ordre de 15 dB, se main- 
tenant encore entre 5 et 10 dB jusque vers 
1000 Hz, et suivie d’une détérioration aux aigués, 
d’autant plus faible que le rayon de courbure est 
plus grand, en ce qui concerne les calottes A, B,C: 
la calotte L se comporte de maniére analogue 
jusque vers 500 Hz, est sensiblement équivalente 
au panneau plan entre 500 et 3000 Hz, pour 
reprendre quelque avantage (de l’ordre de 5 dB) 
au dela. 

Des expériences complémentaires ont montré 
que le comportement restait sensiblement le 
méme, que le son arrivat sur la face convexe ou 
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sur la face concave; et que le serrage réalisé était 
suffisant pour qu’un platrage soigneux du con- 
tour n’apportat pas de modification appréciable 


‘de énergie transmise. 


Ces résultats présentent un certain intérét 
pratique, du fait que ce sont précisément les sons 
graves qu'il est le plus difficile de combattre, et 
qui, une fois les parois de l’enceinte traversées, 
se laissent le plus difficilement absorber par des 
revétements superficiels. D’autre part, leur allure 
est peu habituelle, de sorte qu’il a paru utile d’en 
chercher l’explication théorique. 


2. Théorie approchée des vibrations d’une 
calotte sphérique 


Les équations générales aux vibrations d’une 
enveloppe mince de révolution [2] prennent une 
forme relativement simple lorsqu’il s’agit d’une 
portion de sphére, et que les déformations peu- 
vent étre considérées comme fonctions de la seule 
colatitude §. Elles peuvent alors s’écrire: 

T,+(T,— T,) cot#+ Kk? u= KV, 
KV'+ KV cot@+(T,+T,)+ Kk?w=0, (1) 
G,+(6,—G,) cot6 + KrV=0 
ou les inconnues u (déplacement tangentiel), w 
(déplacement radial), V (cisaillement radial) sont 
liées aux composantes de la tension par les rela- 


tions T, = K[u'+ cucot§—(1+ o)w], 
T, = K[ou'+ ucot§—(1+ a) wu], 


D 
= = [u’+ w"+ o(u+ w’) cot], 


D 
Ga = 7 [o(u'+ w") + (u+ w’) cotO] 


et ot lon a posé 


Kh(1 
E désignant le module de Young, ¢ le coefficient 
de Poisson, m la masse de l’enveloppe par unité 
de surface, et m; une pulsation propre. 
En dérivant une fois la deuxiéme équation (1) 
et en remplacant les tensions par leurs expres- 
sions, les éqs. (1) prennent la forme 


u"+u' cotd—— 5 +(1—o + k?)u=V+(1+o)w’, 

(u+ w’)”+ (u+ w’)'cot@ — wi +(1—o)(u+ w’) 
= tn V, (2) 

6-10 (u's u’ cotd—— + (24+ 20—k?) w’ 
= V"+V'cotd— on 5: 
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Soit L le symbole de l’opérateur différentiel 


linéaire a2 mee, 4 1 
de OO" do sin?” 
et de plus l—o+k?=a,/ 2+4+20—k?=c; 


le systéme (2) prend alors la forme 
(L+ a)u—V—(1+0)w'=0, 


3 
(L+1—o)u— a V+(L+1—o)w’=0, (3) 
(l+o)Lu—LV+cwu'’=0 


et l’on voit qu'il est possible d’y satisfaire en 
prenant u,V,w’ proportionnelles 4 une méme 
quantité et en donnant 4 la lettre L une valeur 
numérique qui annule le déterminant de (3), au 
lieu de la considérer comme représentant un 
opérateur. Autrement dit, on prendra pour V et 
w’ des expressions telles que Au et Bu, et pour 
u une solution de l’équation différentielle 


Lu= Au, (4) 


en choisissant pour la constante 4 une valeur qui 
vérifie ’équation déterminante 


At+a —1 —(1+ 0) 
A~I—o —Kr?/D i4+1—o| =0. (5) 
—(l+o)A —A c 


On reconnait, dans (4), l’équation différentielle 
associée de Legendre du premier ordre ; sa seule 
solution qui soit réguliére au pole 0=0 de la 
calotte est la fonction associée de premiére 


bree u = P}(cos8), 
dont le degré (non entier) » est lié a la valeur de 
la constante / par la relation 
A=—v(v+1). 
I] faut encore satisfaire aux conditions aux 
limites du probléme, lesquelles sont 
u(0)=0, (9) = 05 
pour ce faire, il suffira de prendre pour » l’un des 
nombres de la suite {vi} définie par l’équation 
P},(cos6y) = 0, 
puis de choisir convenablement la constante d’in- 
tégration dans la quadrature qui donne w, soit 
w = B[P, (cos 6) —P, (cos6,)]. 


Reste l’équation déterminante (5), laquelle lie /,, 
choisi comme ci-dessus, aux k; encore arbitraires 
qui figurent dans a et c, et qui est par conséquent 
l’équation aux valeurs propres de la calotte. 
Cette équation peut encore s’écrire 


Kr3/D A+1—oa A+1—o 
A —A(1+¢) c =F 
1 o+a —(1+o) | 
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vu la prédominance numérique du facteur Kr?/D, 
qui est de l’ordre de 4-10° pour les matériaux en 
question ici, les zéros de ce déterminant sont 
sensiblement ceux du mineur relatif a Kr?/D, 
d’out P’équation approchée 
(1+o)?—c(a+/A)=0 
soit encore, puisque o est voisin de 1/3, 
kt + (A; — 2) k? — (0,89 A: + 1,77) = 0. 


On remarquera que cette approximation revient, 
tout en gardant la forme correcte de u,V,w, a 
Mate Ne : : : 
négliger l’énergie de flexion devant celle d’ex- 
tension pour le calcul des fréquences propres 

seulement. 
Pour chaque valeur de /; (donc de vi), ?équa- 
tion bicarrée admet deux racines, sensiblement 


égales a nee) 0,89 vi (vi +1)—1,77 
hae vi (vit 1)+2 
et kay = vi(vit 1) ide 


Le Tableau II indique les premiéres valeurs de 
vi(vi+1) pour différents angles d’ouverture de la 
calotte. 


Tableau II 
A 20° 50° 
Vy 10,75 il Soll 4, 
¥,(%,+1) | 126 56: oh iwuel 20 


Pour les calottes A,B,C, on a K—=7-10!°; m= 
0,27; r=27sin0,); de sorte que les pulsations 
propres sont de l’ordre de 10*k;. Celles du groupe 
II, qui d’ailleurs croissent réguliérement avec le 
rang i, sont toutes dans les aigués, et contri- 
buent peu a la transparence du hublot. Pour le 
groupe I, on constate que les valeurs des k; ten- 
dent a s’accumuler au voisinage d’une borne 
supérieure lorsque 7 croit, et que d’ailleurs cette 
borne n’est point trop élevée; la Fig. 2 permet 
de fixer les idées sur les valeurs numériques, dans 
le cas des calottes étudiées. 

Lorsque la sphéricité de la calotte tend vers 0, 
a restant constant, toutes ces fréquences tendent 


7 02 Odin 


2345 
——= Fréquence 


10 kHz 


Fig. 2. Répartition des fréquences propres d’un hublot 
sphérique en aluminium, de 1 mm d’épaisseur, en 
fonction de l’angle d’ouverture 6). 
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vers 0, et celles qui se manifestent sont désormais 
dues aux vibrations de flexion négligées jus- 
qu’alors, mais que l’annulation de 6), met en 


valeur. Les fonctions propres, de leur cété, ten-' 


dent vers des fonctions de Bessel suivant un 
mécanisme classique, et l’expression, bien connue, 
des pulsations propres du disque plat auquel on 
arrive finalement est 
D @P 
Or = 
ma 
On vérifie que la plus basse d’entre elles est en 
deca de la gamme audible dans le cas envisagé, 
et elles couvrent toute la gamme d’un spectre 
assez serré, puisque l’écartement entre fi et fi +1 
est sensiblement 


Apa" 2 YP fm ayy 


a 2am 


de sorte qu’on en rencontre une tous les 80 Hz 
au voisinage de f= 500 Hz, par exemple. 

Si l’on considére que l’énergie transmise par le 
panneau se présente sous la forme d’une série ot 
chaque terme correspond a une pulsation propre, 
et est affecté d’un diviseur fonction croissante du 
désaccord entre cette pulsation et celle de ’onde 


MEETING OF 1so/TC 43 - ACOUSTICS 


ACUSTICA 
Vol. 4 (1954) 


excitatrice ; que par suite la plus grande partie 
de énergie transmise est celle due aux termes 
quasi-résonants pour lesquels le diviseur est trés 
petit [3]; que de tels termes seront présents a 
toute fréquence dans le cas du disque plan, mais 
n’apparaitront que dans la zone hachurée de la 
Fig. 2 pour une calotte sphérique, on s’explique 
fort bien les résultats expérimentaux donnés plus 
haut. 

La théorie approchée que l’on vient d’établir 
peut également servir de guide dans l’établisse- 
ment de hublots destinés 4 une application parti- 
culiére; mais si l’on. désirait évaluer la trans- 
parence avec quelque précision, il faudrait tenir 
un compte plus exact de la flexion, ce que les 
équations précédentes permettent d’ailleurs de 
faire, au prix de quelques complications supplé- 


mentaires. 
(Regu le 6 avril 1954.) 
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FIRST MEETING OF ISO/TC 43-ACOUSTICS (MUSICAL PITCH) 


On 20th and 21st October 1953 the first meeting of the 
Technical Committee 43 of the International Organization 
for Standardization (ISO/TC 43) was held in London. Its 
most important resolution was the unanimous re-affirma- 
tion of the value of 440 c/s for the Standard Tuning Fre- 
quency, fixed by the 1939 Conference on Musical Pitch. 
The standard tuning frequency is the frequency for the 
note A in the treble stave. 

France agreed to this value of 440 c/s (although there 
were some people, particularly singers, who thought this 
value too high), on the understanding that it would not be 
exceeded and that it would be strictly respected by all 
countries and by all interests. 

The German delegation wished to draw a distinction 
between the standard frequency of 440 c/s and the real pitch 
of performance. They were in favour of requiring all instru- 
ments to be built so that they could be tuned to 440 c/s, 
but they thought no attempt should be made to control 
the pitch used by musicians, 

The Italian delegation said that by law the standard 
frequency in Italy was 435 c/s, but in many cases in prac- 
tice the frequency was higher than this, and was in fact 
nearer to 440 c/s. The ISO should not specify tolerances 
for the pitch during performances, but should only fix the 
frequency for tuning purposes. 

It was suggested that the ISO recommendation should 
be as simple as possible and that any additional informa- 
tion should be given in the national standards. The title 
of the standard was changed to “Standard Tuning Fre- 
quency (Standard Musical Pitch)’. There was conside- 
rable discussion regarding the accuracy to be prescribed 
for tuning devices and finally + 0.5 ¢/s was adopted, 


although instrument manufacturers may find it desirable 
to use tuning devices having a higher accuracy. The French 
delegation proposed a reference to a temperature of 20°C 
but other delegations expressed the view that instruments 
should be capable of being tuned to the standard frequency 
of 440 c/s irrespective of the temperature. 

With regard to the future programme of work there was 
considerable discussion on the question whether all the 
work on acoustics should be dealt with exclusively by 
either the International Organization for Standardization 
(ISO) or the International Electro-technical Commission 
(IEC), or whether the work should be divided between these 
two organizations and carried out with the closest possible 
liaison between the two. 

Finally a resolution was adopted proposing to the 
1SO-Council that ISO/TC 43 should conduct work in 
some purely acoustical fields (1. Acoustical terminology, 
2. Threshold of hearing, 3. Architectural acoustics, including 
a) Sound insulation, and b) Absorption and reverberation, 
4. Airborne noise) and that the meetings of ISO/TC 43 and 
IEC/TC 29 should be held at the same time. 

France, Germany, Italy, the Netherlands, Switzerland, 
the United Kingdom and the United States of America 
had sent delegates, and observers were present from India, 
the European Broadcasting Union, the Technical Committee 
29 of the International Electro-technical Commission (LEC/ 
TC 29, Electro-acoustics) and the General Secretariat of 
the ISO. The Director of the British Standards Institute, 
Mr. BINNEY, welcomed the delegates to British Standards 
House, and Mr. Broapwoop and the B.B.C. entertained 
the delegates in a most delightful manner. 

R. VERMEULEN 


UBER DIE BEWEGUNGSGLEICHUNG UND DEN IMPULSSATZ 


IN VISKOSEN UND KOMPRESSIBLEN MEDIEN 


von F. E. Borenis 


California Institute of Technology, Passadena, Cal., USA. 


Zusammenfassung 


Die Bewegungsgleichung und der hydrodynamische Impulssatz werden in vektorieller Form 
aus dem allgemeinen Ausdruck fiir den Spannungstensor in viskosen Medien unter Einbeziehung 
der Kompressibilitat sowie der Scher- und Volumviskositat hergeleitet. Es folgen einige Bemer- 
kungen tiber die Anwendung einer eindimensionalen Betrachtungsweise auf die durch akustische 
Wellen verursachte Triftstromung in viskosen Medien. 


Summary 

The equation of motion and the hydrodynamic momentum theorem are derived vectorially 
from the general expression for the stress tensor in a viscous medium with due consideration of 
the compressibility and the shear- and bulk-viscosities of the medium. Some remarks are added 
concerning the use of the one-dimensional specialization of the equations for an explanation of 
the hydrodynamic flow caused by acoustic waves. 


Sommaire 

On a établi sous forme vectorielle, 4 partir de l’expression générale du tenseur des efforts, 
léquation du mouvement et le théoréme hydrodynamique de l’impulsion pour les milieux vis- 
queux, en tenant compte de la compressibilité ainsi que de la viscosité de cisaillement et de 
volume. Suivent quelques remarques sur l’emploi d’une méthode unidimensionnelle dans ]’étude 


du courant de déplacement produit dans les milieux visqueux par des ondes acoustiques. 


1. Einleitung 


Die Wechselbeziehungen zwischen akustischen 
Vorgiangen und turbulenter Bewegung in viskosen 
Medien gewinnen in letzter Zeit zunehmend an 
Interesse [1]---[4]. Zur theoretischen Behandlung 
von Problemen dieser Art benétigt man die Euler- 
schen Bewegungsgleichungen in allgemeinerer 
Form als sie gewéhnlich in Lehr- oder Hand- 
biichern zu finden sind..Neben der Scherviskosi- 
tat, die mit der Formanderungsgeschwindigkeit 
eines Massenelementes verkniipft ist, spielt hier 
auch, besonders in Fliissigkeiten, die Volumenvis- 
kositat (second-, dilatational- oder bulk-viscosity) 
eine Rolle, die ein MaB® fiir die viskosen Krafte 
bei einer reinen Volumenanderung eines solchen 
Elementes darstellt. Die durch eine akustische 
Welle in viskosen Medien erzeugte Triftstrémung 
(Quarzwind, hydrodynamic flow) eréffnet eine 
Méglichkeit zur experimentellen Bestimmung des 
Koeffizienten der Volumenviskositat [1], [5]. 

Ein sehr brauchbares Werkzeug bei der Be- 
handlung der genannten Probleme stellt auch der 
hydrodynamische Impulssatz dar, den man ge- 
wohnlich nur fiir nicht viskose Medien hergeleitet 
findet [6], [7]. Uber seine allgemeine Fassung 
unter Beriicksichtigung der viskosen Krafte und 
der Kompressibilitat finden sich in der gebrauch- 
lichen Literatur nur vereinzelte Hinweise. 


Es mag daher nicht iiberflissig sein, eine kurze 
Herleitung der Bewegungsgleichungen und des 
Impulssatzes in allgemeinerer Form zu geben. 
Wir werden uns dabei begniigen, die Viskositat 
als eine Konstante, d. h. von Ort und Dichte un- 
abhangige Gréfe zu betrachten; eine Erweiterung 
fiir variable Viskositat laBt sich unschwer an- 
bringen. 

2. Herleitung der allgemeinen 


Bewegungsgleichung 


Der Spannungstensor in Medien mit innerer 
Reibung besitzt die folgenden Komponenten ([6], 
Si92) 


: Ou 
Txx = —p—A divv + 2u aa? 

: Ov 
Ty = —p—A divy + 2p fal 


Tx = —p—Adivv + ee 
(1) 


T= Tye =p (5 =} =). 
Oy ox 
Ov Ow 
Ty,=T,, erred 
si ls ( Oz ' by 
ow Ou 
re aay ( Ox ‘a Oz ) 
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v (u,v,w) ist der Geschwindigkeitsvektor, p der 
hydrostatische Druck, 4 und y sind Viskositats- 
Koeffizienten. 


Definitionsgemaf folgt daraus der ,,mittlere Nor- 


maldruck*“‘ (2) 
ite : (ee f Tyy+ de) =p—(Fu—2] div v. 


Gewohnlich wird angenommen, das —// mit dem 
hydrostatischen Druck p identisch ist und da} 
dieser nicht von der Expansionsgeschwindigkeit 
divv abhangt. Dies fiihrt zu der bekannten 
StokeEsschen Relation A= - u. Fir einatomige 
Gase ist diese Beziehung korrekt; fiir mehr- 
atomige Gase sind die Abweichungen davon bei 
nicht zu kleinen Drucken wahrscheinlich nicht all- 
zu groB. In Flissigkeiten treten jedoch unter Um- 
standen erhebliche Abweichungen auf. Dies legt 
die Einfithrung eines Koeffizienten der Volumen- 
viskositat yu’ nahe, der proportional der Volumen- 
anderungsgeschwindigkeit eines Massenelementes 
gesetzt wird und im Hinblick auf Gl. (2) tblicher- 
weise durch die Beziehung 


—IT=p—p' divv (3) 


definiert wird!. Damit folgt 


2 4 
A= zpe—p’ oder 2u—A=zutwm’. (A) 


3 3 
Die Eulersche Bewegungsgleichung lautet in 
allgemeinster Form 
Dv ov ° 


0p, =e a te Viv=f+divT, (5) 


1 Das System der Spannungen (1) léBt sich fiir ein- 
atomige Gase aus statistisch-kinetischen Betrachtungen 
herleiten. Unter der Voraussetzung kleiner Abweichungen 
von der Maxwettschen Verteilung ergibt sich fiir solche 


Gase die StokEssche Beziehung 7 = 4 yu, d.h. ein solches 


einatomiges Gas zeigt keine Volumenviskositat (vgl. auch 
die historische Bemerkung von G. G.Sroxes, Math. phys. 
Papers 3 [1851], 136). Die erste Herleitung der Span- 
nungen in einem stark verdiinnten einatomigen Gas mit un- 
gleichmaBiger Dichte- und Temperaturverteilung findet 
sich in einer Arbeit von MAXWELL [10], in der die Gasatome 
als reine StoBzentren mit einem Kraftwirkungsgesetz r-° 
angesehen werden. Von J. E. Jones [11] wurden diese Be- 
trachtungen auf Kugeln mit Eigenvolumen und einem r-”- 
Gesetz erweitert. Die Verhiltnisse in dichten Gasen be- 
trachtet Y. Rocarp [12]. Eine kinetische Theorie fiir mehr- 
atomige Gase miiBte die Existenz weiterer Freiheitsgrade, 
d.h. von Anregungszustinden auf quantenmechanischer 
Grundlage, einbeziehen. Diese Effekte machen das Auftre- 
ten von Relaxationserscheinungen und damit einer an zeit- 
abhiangige Vorgiinge gekniipften Volumenviskositat plau- 
sibel. Eine solche gaskinetisch- quantenmechanische Theorie 
der Viskositat st6Bt auf erhebliche Schwierigkeiten und liegt 
bis heute nicht vor. Uber die phanomenologische Einbezie- 
hung von Relaxationswirkungen vgl. E. Tisza [13] und M. 
Kouter [14]. Eine ausfiihrliche historische Ubersicht findet 
man bei S. CHAPMAN und T.G.Cowtrne [15]. 
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wobei f die Resultierende der an der Volumen- 
einheit angreifenden auBeren Krafte, o die Dichte 
und div T die Tensordivergenz des durch Gl.(1) 
beschriebenen Spannungstensors bedeuten. Durch 


das Symbol div T ist ein Vektor beschrieben, 
dessen x-Komponente gegeben ist durch 


0 Shes. 0 thee a fie 
+" + 


(divT).= On By pe (6) 


Die tibrigen Komponenten erhalt man daraus 
durch zyklische Vertauschung. 
Die Kontinuitatsgleichung lautet 


50h pit, Meas 

agi divov=0. (7) 
Multipliziert man diese Gleichung mit v und 
addiert sie zur linken Seite von Gl. (5), so erhalt 
man 


= (ev) +vdivev+o(v- V)v=£+ div T. (8) 

Der Vektor div T’ laBt sich, wie man durch 
Einsetzen seiner Komponenten aus Gl. (1) in 
Gl. (6) leicht verifiziert, auf die folgende Form 


bringen: 
div T = —grad p+ (uw—A) graddivv+uAy. (9) 


Benutzt man die Vektoridentitat A v= grad divy 
—rotrotv, so kann man auch schreiben: 


3 (10) 

div T= — grad p + (2u—A) grad div v—rotroty. 

Fihrt man Gl. (10) in Gl. (8) ein, so erhalt man 

unter Beriicksichtigung von Gl. (4) die allgemeine 
Euersche Bewegungsgleichung in der Form: 

0 

t 


(ov) +vdivov-+ o(v‘\/)¥ 


in) 


(11) 


4 
= f—gradp+ (50 +W) grad div vy —yrotroty. 


3. Der Impulssatz 


Den Impulssatz gewinnt man am einfachsten 
durch Integration von Gl. (8) tiber ein endliches 
Volumen und Transformation der dabei entste- 
henden Volumenintegrale in Oberflachenintegrale. 


Eine bekannte Verallgemeinerung des GauBschen 
Satzes? liefert 


[eivtav— h(a?) da, 
y 


wobei n die 4uBere Normale zum Flachenelement 


> Siehe z. B. B. FRANK-MieEseEs, Differentialgleichungen 
der Physik. Verlag Vieweg, Braunschweig 1930, Bd. I, 
Kap. 2, §. 4. 
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da bedeutet und (n T) einen Vektor darstellt, 
dessen x-Komponente z.B. gegeben ist durch 


(n sie = Nx Wr! Ny Try + Nz 1 Bs . (12) 


Die iibrigen Komponenten folgen daraus durch 


zyklische Vertauschung. Der Vektor (n T) kann 
nun in folgender Weise dargestellt werden: 


: (13) 
(n T) =—np—nd divv+ 2u(n-V)v+ uw [nx roty]. 


Dies 1a8t sich leicht folgendermafen verifizieren: 
Unter Verwendung der Identitat 


[mn x rotv] = grad (n-v) —(n-\/)v 
erhalt man aus Gl. (13) 
; (14) 
(n T) =—np—nA divv+ u(n:\/)v + wegrad (n-v). 


Bildet man z. B. die x-Komponente von Gl. (14), 
so findet man 


(nT), =—n, p—nx Adivy a 


Dieser Ausdruck ist mit Gl. (12) identisch, wenn 
man dort die Spannungskomponenten des Ten- 
sors in Gl. (1) einftihrt. Analog erhalt man die 
Identitaéten in den beiden iibrigen Koordinaten- 
richtungen. Die linke Seite von Gl. (8) kann nun 
mit Hilfe einer bekannten Beziehung transfor- 
miert werden: 


[lv divov+o(v-V)v] dV= ev(v-n) da. (15) 


Unter Verwendung der Beziehungen (4), (8), (13) 
und (15) erhalt man schlieBlich den Impulssatz 
in der folgenden Fassung: 


© feva V+ bov(v-n) da 


ia 


_ pup da + (uw >") ndivvda ate 


£2 wv )vda+ ud [ax roty] da. (16) 


Die linke Seite stellt die gesamte Impulsande- 
rung innerhalb des von der Integrationsflache 
umschlossenen Gebiets dar; der erste Term die 
lokale, rein zeitliche Anderung des Volumeninte- 
grals iiber die Impulsdichte ov, der zweite Term 
das Oberflachenintegral iiber die Normalkompo- 
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nente des Impulsflusses (Impulsdichte mal Ge- 
schwindigkeit) durch die Oberflache. 

Die rechte Seite enthalt die Resultierende der 
auBeren Krafte F, die Druckwirkung iiber die 
Oberflache und die Oberflachenwirkung der vis- 
kosen Krafte. Die Volumviskositat tritt, wie zu 
erwarten, nur in einem Term w’ divv auf, der die 
Volumenexpansions-Geschwindigkeit enthalt. 

Der Impulssatz in der Form der Gl. (16) bietet 
den Vorteil, dass er im wesentlichen Oberflachen- 
integrale tiber das zu betrachtende Gebiet ent- 
halt. Fir periodische Vorgange verschwindet im 
Zeitmittel die lokale Anderung der Impulsdichte 
und Gl. (16) enthalt nur zeitlich gemittelte Ober- 
flachenintegrale. Die gesamte ,,konvektionelle“ 
Anderung des Impulses im Innern des betrachte- 
ten Volumens V, die mit der raumlichen Ver- 
schiebung von Volumenelementen bzw. Partikeln 
innerhalb V verknipft ist, laBt sich auf ein 
Oberflachenintegral iiber die Normalkomponente 
des Impulsflusses zuriickfiihren. Die gesamte 
Wirkung der viskosen Krafte innerhalb V redu- 
ziert sich auf eine Oberflachenwirkung, die durch 
die entsprechenden Glieder auf der rechten Seite 
von Gl. (16) wiedergegeben wird. Physikalisch ist 
dies durch den Umsiand verstandlich, da die 
viskosen Krafte als innere Krafte sich bei der 
Integration itiber das Volumen paarweise auf- 
heben, so daB in der Summe nur eine Oberflachen- 
wirkung iibrig bleibt. 


4. Bemerkungen zu eindimensionalen 
Betrachtungen 


Wir beschlieBen die allgemeinen Betrachtun- 
gen mit einigen Bemerkungen zum eindimen- 
sionalen Sonderfall, der zwar exakt nicht reali- 
sierbar ist, aber in gewissen Fallen niitzliche 
Schliisse erlaubt. 

Aus der Bewegungsgleichung (11) ergibt sich 
mit 0/dy = 0/éz=0, v=u und bei Abwesenheit 
auBerer Krafte 


(17) 
é rd) ee iG en 
ap (ou) + ae (2¥ bag 2 vale: 
Fir zeitlich periodische Vorgange folgt daraus 


nach zeitlicher Mittelung und Integration nach x 


ou? + p— & p+ u') (= ) =ou®—Tiz—C, (18) 
wobei die Konstante C unabhangig von t und x 
ist. Dieselbe Gl. (18) findet man auch direkt aus 
der eindimensionalen Spezialisierung des Impuls- 
satzes (16) bei zeitlicher Mitteilung und Inte- 
gration iiber den Querschnitt Eins. Bezeichnen 


410 


wir mit pr =p— (5 b+ u') (cu/ox) die totale 
mittlere Druckwirkung in der x-Richtung, so gilt 
im Zeitmittel? 


ou2+ pr = C. (19) 


Nachdem ou? als Dichte des Impulsflusses 
(ou:u) durch den Einheitsquerschnitt interpre- 
tiert werden kann, besagt Gl. (19) physikalisch, 
da jedem Verlust an Impulsflu8 eine Druckzu- 
nahme entspricht. Betrachtet man zwei parallele 
Ebenen x, und xy, so stellt sich nach Gl. (19) bei 
einem ImpulsfluBverlust —A ou? der akustischen 
Welle durch Absorption ein positiver mittlerer 
Druckanstieg Ap; zwischen diesen beiden Ebe- 
nen ein. An der Volumeneinheit halten sich die 
GréBen Ap, und —A ou? das Gleichgewicht. Die 
Absorption der Welle ist in der Naherung erster 
Ordnung durch (= ier uw’) bestimmt, wie aus der 
Bewegungsgleichung (17) hervorgeht; in den 
Druckanstieg gehen daher sowohl die Scher- als 
auch die Volumenviskositat ein. 

Fox und HErzFELpD [9] haben die Beziehung 
(19) zu einer vereinfachten Erklarung der durch 
akustische Wellen in viskosen Medien verursach- 
ten Triftstrémung herangezogen. Ihre Argumen- 
tation geht dahin, da durch die Scher- und 
Volumviskositat der Impulsflu®B eine Verminde- 
rung in der Fortschreitungsrichtung der Welle er- 
fahrt und ein korrespondierender Druckgradient 
die Triftstrémung zustande bringt. Als treibende 
Krafte kann man die Komponenten Tj; in Gl. (1) 
ansehen, die von Scher- und Volumenviskositat ab- 
hangen, wahrend die viskosen Schubspannungen 
Tix, die der Scherviskositaét allein proportional 
sind, der Strémung entgegenwirken. Dies macht 
plausibel, daB die Strémungsgeschwindigkeit, im 


3 Uber eine Anwendung der Gl. (19) auf die Messung 
akustischer Strahlung siehe [8]. 
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Einklang mit Eckarts Ergebnis, vom Verhaltnis 
(5 p+ i) [b= — +w'/u abhangt; ihre Messung er- 


_ laubt daher Aussagen iiber den Wert von w’/u [5]. 


Im streng eindimensionalen Problem, d. h. bei 
unendlichem Strahlquerschnitt, kommt eine Stré- 
mung jedoch nicht zustande, weil sich die treibende 
Kraft —A gu? und der Druckanstieg Ap, gerade 
das Gleichgewicht halten. Dieselbe wird vielmehr 
einer Verminderung von Ap, zuzuschreiben sein, 
die durch die endliche Strahlausdehnung bedingt 
ist. Von wesentlichem EinfluB diirften dabei die 
Geschwindigkeitsgradienten ¢u/éy bzw. 0dv/0x 
sein, die bei einer eindimensionalen Betrachtung, 
bei der rotv von vorneherein verschwindet, gar 


nicht ins Spiel kommen. 
(Eingegangen am 22, September 1953.) 
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A NEW METHOD OF COMPUTING ACOUSTICAL FILTERS 
by W. K. R. Liprertr 


Division of Building Research, Commonwealth Scientific and Industrial Research Organization, Australia 


Summary 

A new method of computing the sound reflection and transmission of acoustical filters which 
does not depend on the image parameter theory is presented. The transmitted and reflected 
sound pressure propagating along a duct with n equal filter sections connected in cascade is 
represented by the characteristic factors of one section which consists of a so-called symmetrical 
discontinuity with a duct piece of arbitrary length connected to the input and output. In loss-free 
cases the sound transmission and reflection can be found with the aid of diagrams. It is shown 
how the new method differs from the image parameter theory and that it is superior in cases 
of a small number of practically loss-free filter sections which are terminated by a non-reflecting 
duct system. The sound reflection and transmission of a filter with a reflecting output duct can 
be easily derived from the fundamental case with a non-reflecting terminal and this method 
can also be applied to filters that are terminated by non-equal sections. 


Sommaire 

On présente une méthode nouyelle pour calculer la réflexion et la transmission du son dans 
les filtres acoustiques, sans avoir recours 4 la théorie du paramétre image. La pression du son 
transmis et réfléchi qui se propage le long d’un conduit comportant n cellules filtrantes identiques 
et montées en cascade, est représentée par les coefficients caractéristiques d’une cellule constituée 
par ce qu’on appelle une discontinuité symétrique, ayant un troncgon de conduit, de longueur 
arbitraire, branché 4 l’entrée et a la sortie. Lorsqu’il n’y a pas de pertes, la transmission et la 
réflexion acoustiques peuvent étre obtenues au moyen de diagrammes. On indique les différences 
existant entre la nouvelle méthode et la théorie du paramétre image; cette nouvelle méthode 
est supérieure dans le cas d’un petit nombre de cellules filtrantes qui sont pratiquement sans 
pertes et qui sont terminées par un systéme de conduit ne présentant pas de réflexions. On peut 
facilement déduire la réflexion et la transmission acoustiques dans le cas d’un filtre avec conduit 
de sortie réfléchissante, du cas fondamental d’une terminaison exempte de réflexions, et cette 
méthode est également applicable aux filtres terminés par des cellules qui ne sont pas identiques. 


Zusammenfassung 

Eine neue Methode zur Berechnung der Schallreflexion und der Schalliibertragung von akusti- 
schen Filtern, die nicht von der Theorie der Wellenparameter abhangt, wird vorgelegt. Der 
Schalldruck des iibertragenen und reflektierten Schalles, der sich in einem Rohr mit n gleichen 
und hintereinandergeschalteten Filterabschnitten ausbreitet, wird durch die charakteristischen 
Faktoren eines Abschnittes dargestellt, der aus einer sogenannten symmetrischen Stofstelle 
mit einem davor- und dahintergeschalteten Rohrstiick beliebiger Lange besteht. In verlust- 
freien Fallen kann die Schalliibertragung und die Schallreflexion mit Hilfe von Diagrammen 
gefunden werden. Es wird gezeigt, wie sich die neue Methode von der Wellenparametertheorie 
unterscheidet und daf} sie im Falle einer geringen Zahl von praktisch verlustfreien Filterabschnitten 
mit einem nichtreflektierenden Rohr als AbschluB iiberlegen ist. Die Schallreflexion und die 
Schalliibertragung eines Filters mit reflektierendem AbschluBrohr kann leicht von dem funda- 
mentalen Fall mit nicht-reflektierendem Abschlu8 abgeleitet werden, und dieses Verfahren kann 
auch auf Filter angewendet werden, die mit ungleichen Filterabschnitten abgeschlossen sind. 


1. Introduction 


The transfer constant of an acoustical filter is 
usually determined with the aid of the so-called 
image parameters [1], [2]. The reason for this is 
that it is easy to find the image transfer constant 
for an arbitrary number of acoustical filter sec- 
tions connected in cascade from that for one 
filter section. The value of the image parameter 
theory, however, is often limited in practice; for 
example, if it is necessary to find the transfer 
constant of an acoustical filter which consists of 
a few sections only and which is terminated at 
the input and output terminals by a constant 


and real impedance instead of the image imped- 
ance. It is true that the effective transfer constant 
can still be approximated quite well by the 
image transfer constant in the middle of the stop- 
bands, as will be shown later, but this is not so at 
other frequencies or in the pass-bands. In such 
cases it is theoretically possible to use correcting 
terms (so-called interaction factors) but this is 
usually impracticable because of the rather com- 
plicated frequency characteristic of the image 
impedance of most acoustical filters. 

The purpose of this paper is to present a new 
method of computing the transfer constant (more 
strictly the characteristic reflection and trans- 
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mission factor) of acoustical filters which does 
not depend on the image parameters and is valid 
in the stop- and pass-bands. A specific advantage 


of this method is that the basic parameters (the ° 


characteristic factors of one filter section) can be 
easily measured, as has been shown in a recent 
paper [3], and that in loss-free cases the fre- 
quency characteristic of the transmission factors 
of n filter sections can generally be determined 
with the aid of diagrams. This simplifies the dis- 
cussion of the effect of the various filter para- 
meters on the frequency characteristics consi- 
derably and makes the new method particularly 
valuable in the practical design of acoustical 
filters. 

A brief study of the relation of the charac- 
teristic factors to the image parameters will show 
that it is easy to represent the image transfer 
constant as well as the image impedance ratio 
by the characteristic factors. This indicates that 
the characteristic factors will be useful in other 
(e. g. electrical) filter problems as well. 


2. The characteristic reflection and transmission 
factors of one filter section 


Fig. 1 shows diagrammatically the structure of 
a symmetrical acoustical filter. The filter consists 
of a duct system of constant cross-section which 
has n equal symmetrical (but otherwise arbitrary) 
discontinuities separated by duct pieces of 
length 1. The term discontinuity is used to de- 
scribe that part of a duct between two parts of 
standard cross-section, in which the area, shape, 
medium, direction or specific acoustic impedance 
of the walls changes. The input and output part 
of the duct are of undefined length and at first 
no reflection is considered from the output duct 
(which has the characteristic impedance Z)). 
Thus the acoustical filter can be divided into n 
equal sections. Each filter section is formed by 
the discontinuity and a duct piece of a length 1/2 
at the input and output. The single symmetrical 
discontinuity has the reference planes 1, 1’ and 
2,2’ and may be specified by the characteristic 
reflection and transmission factor R’ = A’ei~ and 
T’ = B’ei® and the characteristic impedance Z, 
of the input and output duct (cf. [3]). The 
characteristic reflection factor is the complex 
ratio of the sound pressure of the reflected wave 
and the incident wave at the input reference 
plane and the characteristic transmission factor 
is the ratio of the sound pressure of the trans- 
mitted wave at the output reference plane and 
the incident wave at the input plane if the output 
duct is terminated by Z,. The characteristic 
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Fig. 1. Structure of multi-sectional acoustical filters with 
similar sections 


i,i’ = input | reference plane of multi-sectional 
0,0’ = output{ filter (or of one filter section) 
le = input 

Dy Get ont reference plane of S. D. 

S. D. = symmetrical discontinuity. 


factors R, T for one filter section with the 
reference planes i,i’ and o, 0’ (ef. Fig. 1) may 
be called R, = Aei* and T,= Bei and can be 
easily derived from R’, T’ and I, where: 


AL a BER (1) 
p=p'—2nl/, | (2) 


and 
a= a’ —21/A, 


A= wavelength. 


These formulae for filters of different length 
are simple and particularly useful for the design 
if 1 is the parameter. 


3. The characteristic factors for n filter sections 


It will be shown below how the characteristic 
factors for n filter sections can be found from 
R and T. This can be done best by electrical 


analogy. Fig. 2 shows an electrical four-terminal 


i/; lo 
r Electrical ) 
Cat 


circuit analogous to one filter section, sound 
pressure and particle velocity being represented 
by voltage and current respectively!. 


Fig. 2. Electrical four-terminal cir- 
cuit analogous to an acousti- 
cal filter. 


1 It may be noted here that it is possible to represent the 
volume velocity (particle velocity multiplied by cross- 
sectional area of the duct) instead of the particle velocity 
by the electrical current. The resultant changes in the argu- 
ment are trifling because the filter systems which are 
investigated here are nothing but ducts connected in cas- 
cade which have the same cross-section at the input and 
output and at all connecting terminals. Thus a change of 
the analogy as mentioned above would not change the 
characteristic factors but merely the dimension of Z, (or 
of any other impedance) which is irrelevant. 
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For a linear, symmetrical and passive four- 
terminal circuit the following relations between 
the input and output voltages and currents are 


valid 

Uj=4y, +42 21,, Uz=2Zy.1,+2Zy1,- (3) 
The coefficients in eq. (3) form the impedance 
matrix and are related to R and T by 
L443 1—R?+T? Lis ve 


Perper 7, rete 4) 


as has been shown elsewhere [3]. 

If filters have to be connected in cascade, as 
in the present case, it is usual to transform the 
eq. (3) to the following form?: 


U,=aU,—bI,, I,=cU,—al,. (5) 
The coefficients in eq. (5) are sometimes called 


the “transfer” matrix and fulfil the following 
relation [4]: 


a—bc=1. (6) 
It can be easily shown from eqs (3)---(6) that 
© heer? gn (1+ R)?—T? 
2S A sa le 
ye oes (7) 
mas yee, 
2T Lon 
or 
_ __(b/Zo) —¢ Zp 
2a+(b/Z,)+¢Z,’ 


; (8) 


Tr : 
2a+ (b/Z,)+¢Z, 
If two equal symmetrical filters with the coeffi- 
cients a, b, c are connected in cascade, the coeffi- 
cients a,, b,, c. of the resultant filter can be 
found by simple matrix algebra: 


a b a b 


c a c a 


a, b,| | 
x 


or making use of eq. (6) the coefficients are 


b, = 2ab, 


(9) 


|€2 G2 


Oy, = Fxane, “{() 


a,=2a*—1, 


In a similar way it can be seen that the coeffi- 
cients dn, bn, en of n filter sections connected in 
cascade can be represented by the following 


matrix equation: 
2 ( ‘ (11) 
This is equivalent to 
@n =a-a, 1+ V(a2_;—1)(a2—1), (12) 
? The negative sign in eq, (5) is caused by the choice of 


the direction of I, in Fig. 2 which is useful because of 
symmetry with respect to the input and output terminal. 


a 6b 


c a 


lan bn 


| Cn An 
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Bn n Beall 
; =< -//2 a aE ere 12 


Gp 


As a,=a, this yields for two, three and four 
filters connected in cascade 


a, =2a?—1 


a, = 4a3—3a (13) 
a,= 8a*—8a?+1 
ave sa aa a a a 
by C2 
bake Honk ei 
La Sere 
a 4 a?—] (14) 
Ss = 84 _8a8_4a 


If n filter sections are connected in cascade the 
resultant four terminal network is still a sym- 
metrical filter and therefore the coefficients a,, 
bn, Cn are related to the characteristic factors 
Rnr= Ane” and Tn= Bnei’ by a formula iden- 
tical in form to eq. (8). It is 
2g (bn/Zo) —€n Lo 

2an+ (bn/Zo) + enZy ” 

T,=- Z ; 
2 an + (bn/Zo) + ¢n Zo 

With the aid of this formula and eq. (7) and 
(12) the required relation between Rn, T, and 
R, T can be derived. 


This yields for the first four filter sections 
(cf. eq. (13) and (14)) 


R, 
(15) 


1 section: 
R,=R, 
T= Lf; 
2 sections: 
T? 
Re = R (1 ar es] ° 
T2 
i= ee © (16) 
3 sections: 
= T? (1— R? + T?) 
R= R[L+ Rae eT 
T3 
T= (1 — R2)?— R?T? > 
4 sections: 
(1 — R?)? + T?(1—2 R? + T?) 
ioe 2 f 
R,= R f rt [ (1 — R?)3— R?T? (2—2 R2+ a} 2 


Tie 


im (1 — R®)3— R?T?(2—2 R? + T?) © 
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It is obvious that with increasing number n 
simplifications are sometimes possible by repeated 
use of transformations with smaller values of n, 
e.g. the characteristic factors for four sections 
could have been readily derived from those for 


two sections 
i T? pepens 
Ry=R, (14+ 7a) — Ra; 


(17) 


However, even with this in mind, it may seem 
at this stage as if the general formulae, eq. (16), 
are far too complicated and not easy to discuss. 
It will be shown in the next sections that this 
is not true. 


4. Acoustical filters with loss-free discontinuities 


If the filter is loss-free or practically so, as in 
many applications, the general eq. (16) can be 
simplified considerably and can be represented 
graphically. 

It has been shown elsewhere [3] that for a 
loss-free symmetrical discontinuity the following 
relation between the characteristic factors is 
valid: 

A? + B=1, 


pep io ES) 


This is equivalent to the following formula: 


R? = —(1— B?) ei? (19) 


whereby it is possible to eliminate R in the 
characteristic transmission factors T’, of eq. (16). 
As T= Be’, the magnitude Bn as well as the 
phase value f, can be represented by B and f. 
This yields the following set of formulae for loss- 
free cases: 


1 section: 

B, = B, 

B, =B6; 
2 sections: 

B,= = ? 

) B4+ 4(1—B?) cos? 
sin2 
6, = arctan E B+ oe | 


= arctan|— eat tanp |; 
ei 2 cos?B — B? we (20) 
3 sections: 
1 = Seay z ; 
BB cos*p (1 —B)(B?—2 cost) 
r sin3 B+ (1 —B4) sing 
p3 = arctan | oad + 3(1— B2) cosB 
4 cos? B — B? 


= arctan 
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4 sections: 
Bt 
BS + 16 cos? (1 — B?) (B2 —2 cos?) * 
‘ sin4B+2(1—B2)sin2B 
Eig oreten | cos4 B + (1—B?) (4cos 28 + 3—B?) | 
8 cost —4 B* cos? 6 
8 cos#B — 8 B® cos®p + BA tang] 


Cc | 


B, 


= arctan 


This result is significant since there are now only 
two independent variables, and further, the 
values of B, and fn.are periodic functions of B 
and can be represented easily by diagrams, B 
being taken as the parameter. This has been 
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Fig. 3a. B, plotted over f for different values of B. 
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Fig. 3b. B, plotted over f for different values of B. 
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done for 2, 3 and 4 filter-sections in Figs. 3a and 
b, 4a and b, and 5a and b. 

From these figures it can be seen that there 
are minima of B, (indicating stop-bands) at 
values of B= + Nx where N=0, 1, 2, 3, --- and 
B-+1. These minima decrease with a decrease 
in B and with an increase in the number of 
filter sections. 
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Fig. 4a. B; plotted over f for different values of B. 
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Fig. 4b. 8, plotted over f for different values of B. 
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Fig. 5a. B, plotted over § for different values of B. ° 
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Fig. 5b. B, plotted over £ for different values of B. 
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Between two successive stop-bands the value 
of B, increases and reaches the value of unity 
(n—1) times so that there are pass-frequencies 


rather than pass-bands. It can be seen that the 


values of B, as well as the ones of 8, are repeated 
if 8 is increased or reduced by Na except for the 
values of f§, if n is an uneven number. In this 
case fn is periodic with § + 2.Nz as can be seen 
from Fig. 4b. 

If necessary, the values of R, (or An and an) 
can be determined easily from the values of Bn 
and #,, because a relation identical in structure 
to eq. (18) is obviously valid for n filter-sections 
also: 

At + Be=—= 7, 


on —Bn= F2/2. - . (21) 


The easiest way to illustrate the applicability 
of the theory is to refer to actual measurements, 
and therefore a simple measured example will be 


described. 


5. Measured example 


A duct with square internal cross-section of 
S=7.5 em (3 in.) and two rectangular bends 
forming a two-sectional acoustical filter is shown 
in Figs. 6 and 7. One filter section consists of the 
right-angled discontinuity with the input refer- 
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Fig. 6. Duct with two rectangular bends and specification 
of the measured quantities. 


Fig. 7. The measured duct system. 
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ence plane 1, 1’ and the output plane 2, 2’ and of 
a duct piece of a length of 1/2 = 7.5 cm (3 in.) 
connected at the input and output. The whole 
filter has the input and output reference planes 
at i, i’ and o, o’. The duct part with the double 
bends was made of hardwood (1.9 em or %4 in. 
thick) and the input and output duct was made 
of brass (0.95 cm or %% in. thick) with the same 
non-reflecting terminal as described in measure- 
ments at a single right-angled bend [5]. The 
characteristic factors of the double bend can be 
derived from the data of the single bend with the 
aid of the diagrams of Fig. 3a and b without any 
other measurements—this will be described first— 
and as those factors have been measured also it 
will be possible to compare the experimental and 
theoretical results. This provides a good check 
on the theory derived above. 

The characteristic transmission and reflection 
factor T’(=B’e}") and R’ (=A’ei”) of a single 
right-angled bend (with the reference planes at 
1, 1’ and 2, 2’ as shown in Fig. 6) are well known 
from earlier measurements [5] and they can be 
easily transformed for the particular position of 
the reference planes of one filter section in Fig. 6. 
Let the characteristic transmission factor of one 


filter section of Fig. 6 be T= Bei’, then 
B=B', p=B'—2al/a, 
I= 15 em (ef. Fig. 6). 


(22) 


In Fig. 8a and b the values of B, A (= 4A’) 
and of #6’ and f are plotted over the dimension- 
less parameter 6=22S// and as @ is propor- 
tional to the frequency, the particular frequency 
scale is given also. 

The values of B and f at various frequencies 
have been used to find the values of B, and f, 
with the aid of the diagrams of Fig. 3a and b, 
which is a quick and easy method and makes the 
specific value of the diagrams obvious. A, and x, 
have been determined with eq. (21). These com- 
puted values of B,, A,, B. and a, are the dotted 
curves in Fig. 8a and b. It can be seen that the 
resultant frequency characteristics for the various 
factors of the double bend are complex curves 
compared with those of the single bend. 

The characteristic factors of the double bend 
have been measured by using the same technique 
as described in measurements with the single right- 
angled bend. The standing wave in the input duet 
and the progressive wave in the output duct were 
explored with a probe tube microphone and the 
measurements were used to determine the sound 
pressure ratios, distance differences and wave- 
length as is shown in an example in Fig. 6. 
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The magnitude of the characteristic factors for one 
rectangular bend, A and B, and for two bends, A, 
and B, (after Fig. 6), plotted over 6 (and the fre- 
quency f). 

eee By measured 

x x A, measured 

——~ — computed curves. 


COMPUTING ACOUSTICAL FILTERS 417 
t 
02 450_____-900 1350 1800 2250 ch 
Brey xm 
5 2 
‘4 : a, +0 
-05 tt a= r i 0 
Bina) \ 1 
nN A 1X8 At (I) 
[x 2. 
BB, Br = fe yt : & @) 
~1 Apes ea mm I Pe I—-05 
eee 
| ty 4 XB L 
{| . iG “Ee 7) T 
— 2: ee 
als) | ry] La 10 
\ \ 
Ba) | ae | 
xIC %. = x 
-20 te ee et \I -15 
0 02 04 06 08 10x70 


—~ 9 


Fig. 8b. The phase values of the characteristic transmission 
factor of one rectangular bend, f’ and f, for two 
positions of the reference planes and of two bends, 
B. and x, (after Fig. 6), plotted over 6 (and the 
frequency f). 
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Table I 

6 Priies Pr Xe Ba 2, po| %2 — Be 
f hs = 2) rca ie Sa As : B, ot | Ag+ BR) 2 
c/s cm cm cm 
300 | 113 0.133 0.88 01924 |'"—T10 = 5 0.06 —1.35 0.98 Ol, 0.96 —0.58 
600 B12 0.264 0.54 0.72 | — 0.2 28.6 0.30 —1.02 0.94. 1.5] 0.97 +0.49 
800 42.2 0.356 0.31 0.56 16.3 | — 5.2 0.53 —1.46 0.85 —1.98 0.99 +0.52 
1000 34.5 0.438 0.30 0.52 8:3) )|(=—8.5 0.54 —0.04 0.80 —0.51 0.93 +0.47 
1125 31.0 0.487 0.79 0.83 4.5 18.3 0.12 —0.43 0.93 =1.03 0.88 +0.60 
1200 28.1 0.537 0.37 0.52 8.3 8.3 0.46 —1.82 0.76 oe 0.79 —0.50 
1300 26.2 0.576 0.13 0.31 5.4 5.8 0.77 —2.18 0.54 —=1.65 0.88 —0.53 
1400 24.4 0.620 0.06 0.20 3.8 4.4 0.89 —2.38 0.37 —1.84 0.93 —0.54 
1500 22.9 0.660 0.05 0.14 2.4 2.8 0.91 —0.59 0.27 —0.05 0.90 —0.54 
1600 21.2 0.712 0.05 0.14 1.1 1.4 0.91 —0.80 0.27 —0.27 0.90 —0.53 
1700 TOs. 0.770 0.07 0.19 0.3 0.0 0.87 —0.95 0.36 —0.47 0.89 —0.48 
1800 19.0 0.800 0.05 0.12 = = 0.91 = 0.23 — 0.88 a 
In Table I the measured quantities and the Ptr/Pmax = ratio of the transmitted sound 


corresponding characteristic reflection and trans- 
mission factors R,=A,e)*: and T,—B,e)* are 
given for various frequencies and for the input 
and output reference planes at i, i’ and o, o’ of 


Fig. 6. The symbols represent: 


iy, = frequency, 

A = wavelength, 

6 = 2xS/4, 

5 = width of the duct, 

Pmin/Pmax = ratio of minimum to maximum 


sound pressure of the standing 
wave, 


pressure (in the output duct) to 
the maximum sound pressure of 
the standing wave, 


Ay = distance difference between the in- 
put reference plane (i, i’) and the 
first minimum, 

Ax = distance difference between the 


output reference plane (0, oe’) and 
the position in the transmitted 
wave where the phase is z different 
from the phase at the maximum 
between the first and second mini- 
mum. 
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From these quantities the magnitude (B,, A.) 
and the phase values (f,, x) of the characteristic 
factors can be easily found (cf. [3]) and they are 


given in Table I and are plotted also in Fig. 8a . 


and b. It can be seen that these measured values 
agree with the dotted (computed) curves as well 
as can be excepted from the accuracy of the 
measurements. 

The last two columns of Table I show that the 
conditions for a loss-free filter are only approxi- 
mately fulfilled. The approximation is, however, 
good enough to justify the application of the dia- 
grams (which represent the simple eq. (20)) so 
that it is not necessary to use the general eq. (16). 
The difference between the phase values («,—/,) 
of the characteristic factors is either approxima- 
tely —z/2 or +2/2 with the change at frequen- 
cies where the magnitude of the characteristic 
transmission factor (B,) has a maximum. These 
measurements illustrate the working of the pro- 
posed method and its ease of application. 


6. Relation between the characteristic factors and 
the image parameters 


In many practical cases the new method of 
computing acoustical filters is superior to the 
conventional methods and this can be further 
seen by considering the relation between the 
characteristic factors (of one filter section) and 
the corresponding image parameters. Fig. 9 


netsh 


Fig. 9. Specification of the image transfer factor 


Uy Us at Uy 41 
the ace Uy 


(E. C. = electrical circuit). 


shows an infinite number of symmetrical filter 
sections connected in cascade. The impedance 
(image impedance) at the input of each filter 
section is equal to the load impedance at the out- 
put and the ratio U,.,/U, of output and input 
voltage is the same for any filter section. As 
nothing but symmetrical networks are considered, 
it is useful to introduce the impedances Z, and Z, 
of the equivalent lattice type circuit. They are 
related to the impedance matrix eq. (3) by 


Z,=2Zy—LZy., 2g=Zy+ Zp. (23) 
Call the image transfer constant J” or the image 
transfer factor F and the image impedance Z; and 
then the following simple relations exist [4]: 
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pe eae AL Acs » (24) 
U, VZ./Z, ae 

Zi=\/Z,"Z. (25) 


With eq. (4), (23), (24) and (25) it yields 
eer (aaa -T)?—R? 
coth 5 (IT)? FR? 


(1+R)?2—T? 
Li= 2) . 


(26) 


and (27) 
For simplification only the loss-free case will be 


discussed. For one section it follows from (26), 
(27) and (18) or (19) 


YS cosB+ B 
coth > ms Te eRe (28) 
Zi cosa + A + sinB + /1—B? 
= = | ite ——,_ (29) 
Lo cosa—A + sinB —//1— B? 


The significant feature of these formulae is that 
obviously the characteristic transmission factor 
(B and 8) alone is sufficient to represent the 
image transfer constant as well as the image 
impedance ratio (normalized image impedance) 
in a simple way. 

In Fig. 10a and b the magnitude (| F'|) and the 
phase value (gy) of the image transfer factor F 


HES 
=| Bleinaas cosB —B 


y seed at 


has been plotted over { for different values of B. 
There are pass-bands for cos? < B* and stop- 
bands for cos?/>B*. At a cut-off frequency, 
cos? = B?. In a pass-band (cos? < B?) 


(30) 


|F'| =] and y = arccos (“°) (31) 
and in a stop-band (cos? > B?) 
1 B 
|F| = +tanh E arctanh (. sa) | 
cosp } (32) 


and pw a, 


The image transfer factor F,=|Fn\eign for n 
filter sections is easily found from 


| Fn]=(|F|)" and gu=n-p. (33) 


With these formulae and Fig. 10a and b it is 
easy to compare the image transfer factors for n 
filter sections with the corresponding charac- 
teristic transmission factors of Figs. 3, 4 and 5. 
There is one important common feature. In the 
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centre of a stop-band (where /} = + Nz) there is 


very nearly 


[Fal > Bn. (34) 


At other frequencies, however, the difference be- 
tween the image and characteristic factors is 
striking. It is obviously advantageous to use the 
characteristic factors directly instead of the 
image parameters in cases where the acoustical 
filter is terminated by a non-reflecting duct and 
has not too many sections. 
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Fig. 10a. Magnitude (|F|) of the image transfer factor 
plotted over f for different values of B. 
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Fig. 10b. Phase (p~) of the image transfer factor plotted 
over f for different values of B. 


7. Acoustical filter terminated by a reflecting duct 
or by dissimilar filter sections 


The specification of the characteristic factors 
is based on the supposition that the acoustical 
filter is or can be terminated by a non-reflecting 
duct system. From a practical point of view as 
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well as for measuring purposes it is useful to 
consider briefly how the transmission and reflec- 
tion of sound are changed if there is a reflection 
from the output duct (or, as it could be expressed, 
if the filter is terminated by an impedance other 
than the constant and real value of Z,). It will 
be shown that even in such cases the charac- 
teristic factors are helpful for describing the 
changes that occur. 

An easy way for specifying the physical situa- 
tion is given in Fig. 11. In (a) of Fig. 11a symme- 
trical filter is shown diagrammatically with the 
input and output reference planes at i, i’ and 0, o’. 
It is terminated by a non-reflecting duct with the 
characteristic impedance Z,. Such filters can be 
specified completely by the characteristic factors 
R and T and by Z. 


©) ial ‘s 


Non-reflecting 


or terminated by Zo 


Fig. 11. Reflection and transmission factors of a symmetri- 
cal discontinuity with a non-reflecting and a reflect- 
ing output duct (S.D. = symmetrical disconti- 
nuity). 

(a) Symmetrical discontinuity with a non-reflecting 
output duct. 

(b) Reflecting terminal. 

(c) Symmetrical discontinuity with reflecting out- 
put duct. 


In (a) of Fig. 11 R and T are symbolized by 
small arrows which show the direction of sound 
propagation, the sound pressure of the incident 
wave at the input terminal being symbolised by 
an arrow of unit length. In (b) a duct with the 
same cross-section as in (a) of Fig. 11 but with a 
reflecting terminal at s, s’ is shown. The complex 
ratio of reflected to incident sound pressure at s, 
s’ is called o. In (c) of Fig. 11 the filter of (a) is 
terminated by a reflecting duct (as specified in 
(b)) with the length of 1, between the reference 
planes at o, 0’ and s, s’. The values for the reflec- 
tion and transmission factors Ry, Ty and gy as 
shown in (c) of Fig. 11 can be easily found by 
superposing all possible reflected and transmitted 
waves at the different terminals; they are 


oT2e—itthlA 
T—oRe sur 
T 


~ [9 Re stun 


Ru = R+ 
(35) 
Tm 
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(35) 


Ry and Ty are the modified values of R and T,, 


if 90 +0. There is a standing wave in the output 
duct as well as in the input duct and the eq. (35) 
are relatively simple, especially if |o|< 1 which is 
of interest for measuring techniques. It may be 
noted that the reflection factor at the terminal 
s,s’ can be represented equivalently by an imped- 
ance value Z,. It is 
by Ai Ze 
Seely ANity 


(36) 


The formulae of eq. (35) are valid for general 
symmetrical filters and can be used for the 
n-sectional acoustical filters as discussed above 
by taking R, and T, instead of R and T. 

If one compares the sound reflection at s, s’ in 
(b) and at i, i’in(c) of Fig. 11, it is obvious that a 
repeated application of the method of the reflect- 
ing output duct shown in Fig. 11 and eq. (35) 
can be used to describe the sound reflection and 
transmission of filters that are terminated by 
dissimilar filter sections as well. Another syste- 
matic way of modifying the structure of the 
filters of Fig. 1 is to change the reference planes 
(ef. eq. (1) and (2)) at various stages of the trans- 
formation to multi-sectional filters. Therefore, the 
range of acoustical filters covered by the theory 
presented above is much wider than has been 
specified by the structure of Fig. 1. 
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8. Conclusions 


A new method of computing the sound trans- 
mission and reflection at acoustical filters is pre- 
sented. This method does not depend on the 
image parameter theory and is particularly 
valuable for symmetrical filters consisting of a 
small number of filter sections only. In loss-free 
cases the sound reflection or transmission of n 
filter sections can be found with the aid of dia- 
grams from the characteristic factors of one 
section. A measured example shows the practical 
applicability of the-method which is considered 
to be superior to the image parameter method 
in cases where the filter is terminated by a non- 
reflecting duct system. If sound reflection from 
the output duct or a dissimilar filter section has 
to be taken into account, the resultant change 
of the reflection and transmission factors can be 
found in a relatively simple way. 

(Received lst December, 1953.) 
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BOOKREVIEW 


A.E.Greenlees, The amplification and distribution 
of sound; 34 ed. (revised) Chapman and Hall Ltd., 
London 1954, 300 pages, 114 fig., 21 x 15 em?, price 35 sh. 


The book is suitable as a general introduction in the 
field of electro-acoustics. It is easy to read and gives a 
reasonably complete review of the various types of micro- 
phones, pick-ups, amplifiers and loudspeakers. 

However, according to the Introduction, the book should 
contain ‘‘the latest developments in the field of electro- 
acoustics, together with all the items, the worker in this 


field is familiar with”. Taking this aim of the book rigor- 
ously one ought to say that the author has not been suc- 
cessful since the specialist in electro-acoustics does hardly 
find any news in the book. Although the newest develop- 
ments are mentioned, subjects like fine-groove recording, 
magnetic tape recording, loudspeaker columns, ete. are 
hardly dealt with, whereas amplifiers nowadays will have 

negative feedback, more modern tubes ete. 
Nevertheless, the book may be recommended to those 
who want a simple and rapid introduction in a new field. 
J. J. Getux, Hilversum 


AKUSTISCHE MESSUNGEN 


AN EINEM MECHANISCHEN MODELL DES INNENOHRES 


Ein Beitrag zur Hydrodynamik der Cochlea 
von H.-G. DiEstTEL 


III. Physikalisches Institut der Universitat Géttingen 


Zusammenfassung 


Kin mechanisches Modell des Innenohres wird konstruiert, das im Mittel etwa 55fach gegen- 
tiber dem Ohr vergréBert ist. Die Frequenzen werden so transformiert, daB die Geschwindigkeiten 
und die Impedanzen im Modell mit denen des Ohres nahezu iibereinstimmen. An diesem Modell 
werden der Druckverlauf in der Fliissigkeit, die Schnelle der Membran, die Impedanz der Membran 
und die Eingangsimpedanz am ,,ovalen Fenster‘ nach Betrag und Phase gemessen. Die Uber- 
einstimmung zwischen gemessenen und gerechneten Werten ist gut. Die Membranauslenkungen 
und die Fliissigkeitswirbel werden fiir verschiedene Frequenzen photographiert. 


Summary 


A mechanical model of the inner ear with on the average 55-fold exageration of scale has 
been constructed. The frequencies are so transformed that velocities and impedances in model 
and ear are nearly the same. In this model, the pressure in the liquid, the motion and impedance 
of the membrane as well as the entry impedance of the “oval window” are measured in magnitude 
and phase, with good agreement between theory and experiment. The deflections of the membrane 
and the vortices in the liquid at various frequencies are photographed. 


Sommaire 


On a construit, avec un agrandissement moyen de 55/1, une maquette mécanique de l’oreille 
interne. On a modifié les fréquences de facon que les vitesses et impédances soient sensiblement 
les mémes dans la maquette et dans J’oreille naturelle. On a mesuré avec cette maquette l’ampli- 
tude et la phase de Ja variation de pression dans le liquide, de la vitesse et de l’impédance de la 
membrane, et l’impédance d’entrée de la «fenétre ovale». Il y a bon accord entre les valeurs 
mesurées et celles calculées. On a photographié pour différentes fréquences les déviations de la 


membrane et les tourbillons dans le liquide. 
1. Einleitung 


Die Erregung der Hérnerven geschieht durch 
mechanische Schwingungen im Innenohr. Das 
Innenohr besteht, physikalisch gesehen, aus 
einem. von starren Wanden umgebenen mit 
Flissigkeit gefiillten Kanal, der Schnecke. Er ist 
durch eine elastische Trennwand, den Ductus 
cochlearis, seiner Lange nach in zwei Teile, die 
Seala vestibuli und die Scala tympani, geteilt. 
Der Ductus cochlearis schwingt als einheitliches 
Ganzes und besteht im wesentlichen aus der 
Basilarmembran, der ReiBnerschen Membran 
und der dazwischenbefindlichen Endolymphe. 
Der Kiirze halber soll in Zukunft mit dem Aus- 
druck Basilarmembran der ganze Ductus coch- 
learis gemeint sein. 


Am hinteren Ende stehen die beiden Scalen 
durch eine Offnung, das Helicotrema, mitein- 
ander in Verbindung. Am vorderen Ende werden 
sie durch je eine elastische Membran, das ovale 
bzw. runde Fenster, abgeschlossen. Wird durch 
den Mechanismus Hammer-AmboB-Steigbiigel am 


ovalen Fenster das Innenohr zu Schwingungen 
angeregt, so pflanzen sich Druckwellen in der Scala 
vestibuli fort. Sie treten iiber die Basilarmem- 
bran in die Scala tympani ein und laufen zum 
runden Fenster. Bei sehr tiefen Frequenzen findet 
der Druckausgleich tiber das Helicotrema statt. 

Es ist v. BEKEsy gelungen, das dynamische 
und statische Verhalten des Innenohres ein- 
gehend zu untersuchen [1]. Er beobachtete langs 
der Basilarmembran fortschreitende Wellen, die 
an einer von der erregenden Frequenz abhangigen 
Stelle ein Amplitudenmaximum aufweisen und 
dahinter schnell abklingen. Ihre Wellenlainge 
wird mit der Entfernung von den Fenstern stetig 
kleiner. Weiter stellte v. BEKEsy fest, daB die 
mechanischen Eigenschaften der Basilarmembran 
stark ortsabhangig sind. Die Federung der Mem- 
bran nimmt mit der Entfernung von den Fenstern 
zu. Das logarithmische Dekrement der Membran- 
schwingungen ist tiber die ganze Membran nahe- 
zu konstant etwa 1,6. Die Dampfung ergibt sich 
also beinahe aperiodisch und erklart die kurze 
Einschwingzeit des Ohres. 
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Die experimentellen Ergebnisse von v. BEKESY 
sind von verschiedenen Autoren theoretisch ge- 
deutet worden. In einer grundlegenden Arbeit 
gab RANKE [2] eine Theorie der Hydrodynamik 
der Schneckenfliissigkeit. Er behandelte die 
Potentialschwingungen als Randwertproblem. 
Die Konstanten der Membranimpedanz wurden 
den Verhaltnissen im Ohr angepabt. 

ZwisLocki [3], PETERSEN und Bocerr [4] und 
FLETCHER [5] beschrankten sich auf den Fall, 
daB die Wellenlange im Kanal grof gegen die 
Querabmessungen ist. Sie erhalten im wesent- 
lichen die gleiche Differentialgleichung, interpre- 
tieren sie aber in verschiedener Weise. Bestim- 
mend fiir die Wellenausbreitung im Innenohr ist 
stets die stark ortsabhangige Federung der 
Basilarmembran. 

Die charakteristischen Membranschwingungen 
wurden von v. BEKESy und anderen Autoren 
auch in Modellversuchen beobachtet [1], [3], [6]. 
Der Druckverlauf im Ohr ist bisher nicht gemes- 
sen worden. BoceErt [7] bildete das Innenohr mit 
Hilfe einer elektrischen Ersatzschaltung nach, 
um den Einflu8 der Dampfung zu untersuchen. 

Das Ziel der vorliegenden Arbeit ist, ein stark 
vergréBertes Ohrmodell zu konstruieren, das den 
Verhaltnissen im Innenohr qualitativ entspricht. 
An diesem Modell sollen die interessierenden 
GréBen: Druckverlauf in der Flissigkeit, Phasen- 
geschwindigkeiten, Auslenkung der Membran, die 
Impedanz der Membran und die Eingangsimpe- 
danz des Modelles gemessen und mit der Theorie 
verglichen werden. 


2. Zur Hydrodynamik des Innenohres 


Im folgenden wird eine kurze Abhandlung der 
Hydrodynamik des Innenohres in Anlehnung an 
die Arbeiten von RANKE, ZWISLOCKI, PETERSEN 
und Bocerr, und FLETCHER gegeben. In der 
Nahe der Fenster ist die Wellenlange groB gegen 
die Querabmessungen der Scalen. In diesem Ge- 
biet wird die Wellenausbreitung durch eine Diffe- 
rentialgleichung beschrieben, die fiir den Fall 
kleiner Amplituden abgeleitet werden soll. FLET- 
CHER beriicksichtigt die Zahigkeitsverluste in der 
Fliissigkeit, kann sie jedoch in der Rechnung ver- 
nachlassigen, da sich ihr EinfluB als gering her- 
ausstellt. Er unterteilt die Basilarmembran in 
kleine quadratische Membranelemente, die gleich- 
phasig und ohne seitliche Kopplung schwingen 
sollen. Die mitschwingende Fliissigkeitsmasse 
wird zur Membranmasse hinzuaddiert, sie ist als 
starr mit dem Membranelement verbunden zu 
denken. Wenn p den Druck in der Fliissigkeit be- 
deutet, wirkt auf ein Element die Kraft pb?. (Die 
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Breite des Elements sei 6.) Die Impedanz eines 
Membranelements, d. h. das komplexe Verhaltnis 
von Kraft zu Schnelle senkrecht zur Membran 


“ ]aBt sich durch den Ausdruck 


: 1 

Z(x)=joM+R+ pew 

der urspriinglich nur fiir punktformigen Kraft- 

ansatz definiert ist, beschreiben. M ist die Ge- 

samtmasse, R der Reibungswiderstand und F die 
Federung. 

Es ist fiir die Messungen einfacher, mit einem 
Einkanalmodell zu arbeiten. Man denke sich die 
Flissigkeit aus der Scala tympani entfernt und 
hat dann eine Fliissigkeitssaule, die an drei Seiten 
schallhart begrenzt ist, wahrend an der vierten 
Seite die Randbedingung durch den Widerstand 
der Membran gegeben wird. An den physikali- 
schen Verhaltnissen andert sich durch den Uber- 
gang zum Einkanalmodell nichts Wesentliches. 
Deshalb wurde diese Form fiir Modellversuche 
meist verwendet [1], [3]. 


Anregung Membran _Helicotrema" 


i Xz t 
Mitty Ly 


Abb. 1. Das Einkanalmodell (schematisch). 


Die Abb. 1 zeigt das Einkanalmodell in sche- 
matischer Darstellung. Betrachtet man einen 
Streifen von der Breite dx, so gilt fiir ihn nach 
dem Newtonschen Grundgesetz bei Vernach- 
lassigung der Zahigkeit: 


jogu——7—. (1) 


Es bedeuten: 

u,v = x- baw. y-Komponente der Fliissigkeits- 
geschwindigkeit, 

P = Schalldruck, 

o@ = Dichte, 

S = Querschnitt des Kanals, 

b = Breite der Membran, 

h =  Tiefe des Kanals, 


wo = Kreisfrequenz. 


Der Druck ist im Falle groBer Wellenlangen iiber 
den Querschnitt konstant. Die mitschwingende 
Fliissigkeitsmasse ist klein gegen die unter einem 
Membranelement befindliche Fliissigkeitsmasse. 
Wegen der Nachgiebigkeit der Membran ist es er- 
laubt, die Fliissigkeit als inkompressibel zu be- 
handeln. Dann ergibt sich fiir die Kontinuitats- 
gleichung 
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Ox 


v ist die mittlere Geschwindigkeit der Membran 
in y-Richtung. Die Impedanz der Membran liefert 
die noch fehlende Beziehung zwischen Druck und 


Schnelle. Es ist 


+5b=0. (2) 


pen 

poe (3) 
Der Faktor n beriicksichtigt die seitliche Ein- 
spannung der Membran und hat nach FLETCHER 
Werte von 0,3---1. Aus den Gl. (1)---(3) folgt die 
Differentialgleichung 


@p 1 0S cp _ jweb'n 
eee Ot OX ZS 


Sie entspricht der Trichterwellengleichung [9], 
[10], jedoch ist hier der Faktor von p eine Funk- 
tion des Ortes. Wegen des komplexen Verlaufes 
von Z wird es im allgemeinen nicht méglich sein, 
die Gleichung in geschlossener Form zu lésen. 
Sind Z, S und b als Funktionen von x bekannt, 
so ist eine numerische Integration méglich. In 
der Nahe einer Stelle x) kann der Querschnitt 
durch den Ansatz S=S,+ 2a(x—x,) angenahert 
werden. In diesem Bereich gilt dann die Glei- 
chung 


o— 


p=0.. (4) 


Pp 2a dp jweben 

0x? Jey Peale As) 
in welcher die Koeffizienten als konstant ange- 
sehen werden kénnen. 

Da im Ohr fiir die meisten Frequenzen keine 
Reflexion am Helicotrema eintritt, soll auch im 
Modell von Reflexionen abgesehen werden. Der 
Ansatz p = p,exp(—yx) liefert nach Einsetzen in 
Gl. (5) fiir die Fortpflanzungskonstante y = x + jf 
den Wert 


wig a\? jwoebin 
r=s+V(s)+*s"- 


Die Fortpflanzungskonstante ist eine komplexe 
GréBe. Man bezeichnet ihren Realteil « als 
Dampfungskonstante, den Imaginarteil / als 
Phasenkonstante. Die Fortpflanzungskonstante 
hangt wesentlich von der Membranimpedanz ab. 
Daneben enthilt sie einen reellen Anteil a/S, der 
von der Querschnittsinderung herriihrt. Da im 
Ohr der Einflu8 der Querschnittsinderung gegen- 
iiber der Anderung der Membranimpedanz nur 
eine geringere Rolle spielt, soll bei der Diskussion 
einiger Sonderfalle einfachheitshalber der Quer- 
schnitt als von x unabhingig angenommen wer- 
den. 


p=9, (5) 


a) Die Membranimpedanz sei ein reiner Fede- 


rungsWwiderstand. Manerhalt y=j | oF b'n/S=jp. 
Die Kopplungsmassen odx (vgl. spater das Er- 
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satzschaltbild Abb. 2) stellen mit der Federung 
Fb’ndx/S ein schwingungsfahiges Gebilde dar, 
in welchem sich ungedampfte, fortschreitende 
Wellen ausbreiten. Ihre Phasengeschwindigkeit 
sitzt die Membranimpedanz zusatzlich eine reelle 
Komponente, so werden die fortschreitenden 
Wellen gedampft. Die Dampfungskonstante bleibt 
stets kleiner als die Phasenkonstante. 


pax ax 
eeule? A , 


U ytdx 
=e 
Py fi b°ndx Cr OnY 
S 
R — 
Vbdx bendx 
Seat 
o—— oO 


Abb. 2. Ersatzschaltbild eines Elementarstreifens. 


b) Die Membranimpedanz sei reell. In 
_14j y/ oot 
Lares 


werden Dampfungskonstante und Phasenkon- 
stante gleich grof. 


c) Die Membranimpedanz sei ein reiner Massen- 
widerstand. Es folgt y=|ob3n/MS=«. Die 
Phasengeschwindigkeit ist unendlich groB. Sie 
bleibt endlich, wenn die Impedanz eine zusitz- 
liche reelle Komponente hat. Die Phasenkon- 
stante ist stets kleiner als die Dampfungs- 
konstante. 

Fir das Einkanalmodell mit ortsunabhangigem 
Querschnitt laBt sich ein elektrisches Ersatz- 
schaltbild des Elementarstreifens der Breite dx 


angeben (Abb. 2). 


Es entsprechen: 


Druck p der Spannung U 
Schnelle u dem Strom I 

Dichte 0 der Induktivitaét/em L 
mechanische der komplexen 


Ableitung/em b?n/ZS Ableitung/em G 


In der Leitungstheorie erhalt man fiir die Fort- 
pflanzungskonstante y= VRG, wobei R den 
komplexen Langswiderstand und G die Ableitung 
bedeuten. Diese Formel erweist sich als mit 
Gl. (6) identisch, da wegen der Vernachlassigung 
der Zahigkeit R=jwL“jwo zu setzen ist. Die 
Beriicksichtigung der Zahigkeitsverluste hatte 
durch Einschalten eines ohmschen Widerstandes 
in Serie mit den Langsinduktivitaten (0/2)dx zu 
erfolgen. 
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Das Ersatzschaltbild erméglicht es, die Schwin- 
gungsverhaltnisse einfach zu ibersehen. Man 
denke sich eine groBe Anzahl solcher Vierpole, 


die sich in der Dimensionierung der Querverbin- ' 


dungen unterscheiden, hintereinandergeschaltet 
und von einer Spannungsquelle gespeist. In der 
Nahe der Anregungsstelle sei die Resonanzfre- 
quenz der Querzweige hoch und werde mit wach- 
sender Entfernung kleiner. Eine Spannungswelle 
von gegebener Frequenz lauft etwa bis zur Reso- 
nanzstelle in die Kette hinein. Dabei steigt ihre 
Dampfung wegen des mit wachsender Entfernung 
geringer werdenden kapazitiven Anteils der 
Membranimpedanz an, wahrend die Phasenge- 
schwindigkeit sinkt. Der durch die Querwider- 
stande flieBende Strom, welcher der Membran- 
schnelle entspricht, ist einerseits durch die an- 
liegende Spannung, andererseits durch die GréBe 
des Widerstandes bestimmt. Obwohl die Span- 
nungswelle gedimpft ist, steigen die Stréme zu- 
nachst an, da die Betrage der Querwiderstande 
abnehmen. An einer von der Frequenz abhangi- 
gen Stelle flieBt ein maximaler Querstrom, wah- 
rend hinter dieser Stelle wegen der Dampfung 
und des Anwachsens der Querwiderstande die 
Stréme sehr klein werden. Die Dampfung der 
Spannungswelle bewirkt, daB die Stelle maxima- 
len Stromes ein wenig von der Resonanzstelle in 
Richtung zur Anregungsstelle hin verschoben ist. 

Fiir den mechanischen Eingangswiderstand des 
Kanals p/u erhalt man analog zu der bekannten 
Formel aus der Leitungstheorie Z..=/R/G@ den 
Ausdruck 


Zmech =VjooZS/b3n. (7) 


In der Nahe der Anregungsstelle ist die Mem- 
branimpedanz fiir tiefe Frequenzen nahezu eine 
reine Federung. Fiir diese Frequenzen ist der 
Eingangswiderstand reell. 

In der Nahe der Resonanzstelle der Membran 
kommen die Wellenlangen in die GréBenordnung 
der Querdimensionen. Die bei der Ableitung der 
Differentialgleichung (4) gemachteVoraussetzung, 
daB der Druck tiber den Querschnitt konstant ist, 
ist nicht erfiillt. Die mathematische Behandlung 
des Problems bietet groBe Schwierigkeiten. 
RANKE hat unter’ der Voraussetzung, dafB die 
Feldgré8en nur von zwei Dimensionen abhangen, 
das Geschwindigkeitspotential in der Form 
F(x + jy) angesetzt und die Lésung fiir den Fall 
groBer Kanaltiefe diskutiert. Die Behandlung 
des Problems hat Ahnlichkeit mit der mathe- 
matischen Behandlung der Oberflachenwellen 
von Fliissigkeiten [8]. 

Gegeben sei ein rechteckiger Kanal der Tiefe 
y=h=const, der an drei Seiten schallhart be- 
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grenzt ist. Die Randbedingung an der vierten 
Seite wird durch die nur von x abhingige Mem- 
branimpedanz bestimmt. Es gelten bei Vernach- 
lassigung der Zahigkeit die Bewegungsgleichun- 
gen 

joov=—épléy, (8) 
ferner die Inkompressibilitatsbedingung 

du/0x + dv/dey = 0. 

Durch Zusammenfassen dieser Gleichungen er- 
halt man A p= 0. Von einer Reflexion am Kanal- 


ende sei wieder abgesehen. Der Ansatz p(x,y) = 
P(y) exp(—y x) fihrt auf 


P(%,¥) = Po 608 (yy) exp (—y%). (9) 
Das Sinusglied in der allgemeinen Lésung der 
y-Abhangigkeit fallt wegen der an der schall- 
harten Wand erfolgenden Spiegelung fort. Die 
komplexe Ausbreitungskonstante y bestimmt sich 
aus der Randbedingung an der Membran. Das 
Verhaltnis Druck zu wandnormaler Schnelle sei 
als Membranimpedanz definiert. Von seitlicher 
Kopplung soll abgesehen werden. Aus den Gl. (9) 
und (8) ergibt sich 


a= (2) ive 
Oo). yeh) 

Nach Erweitern mit h erhalt man die transzen- 

dente Gleichung fiir die Fortpflanzungskonstante 


joou=—éop/ox, 


(10) 


_ jwgh 

yh te rh) ae 

Sie ist aus der Theorie der Schallausbreitung in 
einem mit Schluckstoffen ausgekleideten recht- 
eckigen Kanal bekannt [10], [11]. Morse hat 
konforme Abbildungen angegeben, die es ermég- 
lichen, zu einer vorgegebenen Impedanz die Fort- 
pflanzungskonstanten zu ermitteln. Es gibt stets 
eine unendliche Schar von Lésungen, die benétigt 
werden, um einen beliebig vorgegebenen Anfangs- 
zustand beschreiben zu kénnen, jedoch inter- 
essieren nur die niedrigsten Werte, da die héhe- 
ren starker gedampften Wellenziigen entsprechen. 
In einer geniigend kleinen Umgebung einer 
Stelle x) kann die Impedanz als konstant ange- 
sehen werden. Dann ist die Fortpflanzungskon- 
stante aus Gl. (10) berechenbar. Nach Gl. (9) ist 
der Druck auch von y abhingig; die Ebenen kon- 
stanten Druckes sind gegen die y-Achse geneigt. 
Fiir den Fall groBer Wellenlangen kann in 
Gl. (10) der Tangens durch sein Argument ersetzt 
werden. Dann geht die Fortpflanzungskonstante 
in die des eindimensionalen Problems iiber, wenn 
man beriicksichtigt, das die Impedanzen in bei- 
den Fallen etwas verschieden definiert wurden. 
Sie unterscheiden sich um die Dimension einer 


Flache. 


(10) 
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3. Ahnlichkeitskriterien 


Fiir Modellversuche ist es haufig zweckmaBig, 
die Abmessungen der urspriinglich betrachteten 
Anordnung in einem bestimmten Verhiltnis zu 
verandern. Ebenso wird man bestrebt sein, die 
Frequenzen in ein meBtechnisch bequem zugang- 
liches Gebiet zu legen. 

Die beiden Anordnungen sollen sich in ent- 
sprechenden GréfSen nur in den beiden Mafein- 
heiten unterscheiden. 

Also L=kl,, fi=Cfz, m=9%2> 01=€02 

(l = charakteristische Lange, f = Frequenz, 7 

= Zahigkeitskoeffizient, 9 = Dichte). 

Sollen in beiden Systemen die gleichen Stré- 
mungs- und Druckgesetze gelten, so sind nach der 
Reynoupschen Ahnlichkeitstheorie [8] folgende 
Beziehungen zu erfiillen: 


Plt 
olif? Oo 13 fP ee 
und 
4 Ne 
= . 12 
ah It Oofo ls (12) 


Gleichbedeutend mit Gl. (11) sind folgende Trans- 
formationen fiir die Komponenten M, R und F 
der Membranimpedanz: 


M, M, 


rl ; 13 

ere (13) 
as 

= ; 14 

OfR .. ah aa 

F,0,ftl? = F,0,f212. (15) 


Die Dichteinderung der Fliissigkeit geht nur in 
der ersten Potenz ein und ist daher von geringem 
EinfluB. Bei vorgegebenem VergréBerungsver- 
haltnis « wird die Masse durch Gl. (13) bestimmt. 
Durch geeignete Bemessung der Federung F 
kann die Gl. (15) fiir eine beliebige Frequenz- 
transformation ¢ erfillt werden. Der Reibungs- 
widerstand R ist dann so zu dimensionieren, da} 
die Gl. (14) befriedigt wird. 

Die Beriicksichtigung der Zahigkeit in der 
Perilymphe verlangt zusatzlich die Erfillung der 
Gl. (12). Das kann durch Wahl einer Fliissigkeit 
mit geeignetem Zahigkeitskoeffizienten gesche- 
hen. Wie von FLETCHER nachgewiesen wurde, ist 
die Zahigkeit in der Perilymphe vernachlassigbar. 
Die Erfillung der Gl. (12) ist daher nicht kritisch, 
solange auch bei den Fliissigkeiten des Modelles 
die Zahigkeit vernachlassigt werden kann. 


4, Aufbau der Apparatur 


Das Ohrmodell ist als Zweikanalmodell ausge- 
fiihrt, in den unten beschriebenen Messungen 
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wurde es jedoch als Einkanalmodell verwendet. 
Das Modell besteht aus einem Plexiglasgehause 
von 1 em Wandstarke. Die Grundflache ist trapez- 
férmig, die Breite nimmt mit der Entfernung von 
der Anregungsstelle zu (vgl. Abb. 3, 4 und 14). 
Die auBeren Abmessungen des Modells betragen 
38x (7—12)x9 cm’. Die Stirnwand wird durch 
zwei Gummimembranen yon 1 mm Starke ge- 
bildet, analog den beiden Fenstern des Innen- 
ohres. In Abb. 3 erkennt man links die zur An- 
regung des Modells dienende Stange, die von 
einem elektrodynamischen System  getrieben 
wird und an der Anregungsmembran befestigt ist. 
An der Riickwand befindet sich der Einfillstut- 
zen fiir die Fliissigkeit und die Gummidichtung 
fiir das Mikrophon. Das Gehause kann mit einem 
Deckel verschlossen werden. 


Abb. 3. Das Ohrmodell mit magnetischem Abnahmesystem 
und Wasserschallmikrophon; seitliche Ansicht. 


Die Membran verlauft in der Mitte des Modells. 
Sie besteht aus einer 1 mm starken Gummiplatte 
und liegt ungespannt zwischen zwei Metallein- 
satzen von je 5 mm Starke, die durch Schrauben 
leicht gegeneinander gedriickt werden. Die Ein- 
sitze haben eine Ausfraisung, durch welche die 
Form der Membran bestimmt ist. Sie werden 
durch Plexiglasstreifen gehalten, die von auBen 
durch Schrauben angezogen werden kénnen. 
Diese Schrauben sind in den Abb. 3, 4 und 16 
(hier als schwarze Flecken) zu sehen. Die Mem- 
bran ist trapezférmig. Ihre Linge betragt 32 cm, 
die Breite in der Nahe der Anregungsstelle 0,72 cm 
und am ,,Helicotrema‘‘ 3 cm. In Abb. 4 ist die 
Membran als dunkle Flache zu sehen. Sie tragt 
langs ihrer Mittellinie eine Reihe weiBer Scheib- 
chen. 

Zur Fiillung des Modells wird entgastes Wasser 
verwendet. Es ist darauf zu achten, dal} keine 
Luftblasen im Modell zuriickbleiben, da sie Fede- 
rungsnebenschliisse zur Membranimpedanz dar- 
stellen. Die Fiillung des Modells mit einer Glyze- 
rin-Wasser-Mischung brachte keine Anderung 
gegeniiber der Wasserfiillung. 
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Zur Messung der Schnelle der Membran wur- 
den auf ihr in 1 cm Abstand voneinander kleine 
Plattchen aus einer Eisen-Nickel-Legierung an- 


gebracht. Diese stéren die Membranbewegung ° 


nicht und erzeugen bei ihrer Bewegung durch 
Induktion eine Spannung in einem Abnahme- 
system. Dieses besteht aus einem kleinen Perma- 
nentmagneten von 5mm @, der eine Wicklung 
von 300 Windungen Kupferdraht tragt. Der Ma- 
gnet ist durch eine Schraube in seinem Abstand 
relativ zur Membran verstellbar. Die horizontale 
Verschiebung geschieht mit Hilfe einer Schwal- 
benschwanzfihrung, die in dem geschlitzten Dek- 
kel des Modells lauft. In Abb. 3 erkennt man auf 
dem Deckel das Tiirmchen mit der Verstell- 
schraube, an dem der Magnet befestigt ist. Die 
Empfindlichkeit der Anordnung ist so groB, dab 
der Abstand zwischen Magnet und Membran 
4-5 mm betragen kann. Bei diesem Abstand 
konnte keine Riickwirkung auf die Membranbe- 
wegung festgestellt werden. 


Abb. 4. Das 


Ohrmodell 
Deckel abgenommen; im Vordergrund Mikrophon- 
einfiihrung und Einfiillstutzen. 


vom ,,Helicotrema‘‘ gesehen. 


Die Messung des Druckes geschieht mit einem 
kleinen Turmalinkristallmikrophon von 7 mm @. 
Es ist an einem diinnen Rohr befestigt und von 
auBen verschiebbar. Eine Gummidichtung an der 
Riickwand verhindert das Austreten von Fliissig- 
keit und ein Ubersprechen von Kérperschall auf 
die Stange des Mikrophons. Da der Druck beim 
Einkanalmodell am _ ,,Helicotrema‘* zusammen- 
bricht, bewirkt diese Gummidichtung keine Ver- 
falschung der Verhaltnisse. 
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Das elektrodynamische Anregungssystem wird 
von einem Tieftongenerator iiber einen Verstar- 
ker gespeist. 

Empfangsseitig werden ein Vorverstarker, 1/,- 
Oktavfilter und ein Anzeigeverstarker verwendet. 
Die angezeigte Spannung wird aufberdem mit 
einem Neumann-Pegelschreiber in eine gleich- 
phasige Spannung konstanter Amplitude ver- 
wandelt, die an die Platten einer Braunschen 
Réhre gelegt wird. An das andere Plattenpaar 
schaltet man eine Spannung, die vom Tiefton- 
generator geliefert wird und mit einem DEsGuIN- 
schen Phasenschieber in bekannter Weise phasen- 
verschoben werden kann. Durch Einstellen der 
Lissajousfigur auf Phasengleichheit kann die 
Phasenlage der am Modell gemessenen Spannung 
bestimmt werden. 


5. MeBergebnisse 
a) Druckverlauf in der Fliissigkeit 


Der Druckverlauf wurde nach Betrag und 
Phase relativ zum Membrananfang (x= 1 ecm) 
gemessen. Er ist in den Abb. 5 und 6 eingetragen. 
Die Abszisse stellt die Entfernung von der An- 
regungsstelle dar; Parameter sind die Frequenzen. 

Die Druckamplituden nehmen etwa linear mit 
dem Abstand von der Anregungsstelle ab und 
sind an der Resonanzstelle der Membran im 
Mittel auf 20% ihres Anfangswertes abgefallen. 
Eine Unstetigkeit in der Nahe der Resonanzstelle 
wird nicht beobachtet. 

Die in den Abb. 5 und 6 eingezeichneten punk- 
tierten Linien sind rechnerisch ermittelte Werte. 
Sie wurden mit den in noch zu besprechender 
Weise ermittelten Konstanten der Membranim- 
pedanz durch numerische Integration der Diffe- 
rentialgleichung (4) berechnet. Diese Auswertung 
wurde fiir die drei reprisentativen Frequenzen 
32, 64 und 160 Hz vorgenommen. Die Uberein- 
stimmung zwischen Theorie und Experiment ist 
gut. Mit gleicher Genauigkeit hatten auch in den 
folgenden Abbildungen die berechneten Werte 
eingezeichnet werden kénnen, doch wurde davon 
aus Griinden der Ubersichtlichkeit abgesehen. 

Die Phasenkurven des Druckes in den Abb. 5 
und 6 entsprechen fortschreitenden Wellen, deren 
Phasengeschwindigkeit mit wachsender Entfer- 
nung sinkt. Da das logarithmische Dekrement 
langs der Membran nahezu konstant ist, und die 
Querschnittserweiterung des Kanals nur einen 
geringen EinfluB auf die Fortpflanzungskonstante 
hat, erreicht jeder Wellenzug an seiner Resonanz- 
stelle den gleichen Wert fiir die Phasenkonstante. 
Die Wellenlange der Druckwelle ist danach fiir 
alle Frequenzen an der Resonanzstelle etwa 18cm. 
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Die Abb. 7 zeigt den Verlauf der Phasengeschwin- 
digkeit als Funktion der Entfernung von der An- 
regungsstelle. Die Werte betragen etwa 45 m/s 
in der Nahe des Membrananfangs und sinken bis 
auf wenige m/s in der Nahe des Membranendes. 
Man erkennt, dafs praktisch keine Dispersion auf- 
tritt. 

Aus dem gemessenen Druckabfall laBt sich die 
Dampfungskonstante « ermitteln. Sie ist in Abb.8 
in dB/em als Funktion der Entfernung und der 
Frequenzen eingezeichnet. Sie erreicht fiir alle 
Frequenzen an der Resonanzstelle etwa den glei- 
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chen Wert x«=3dB/cm, entsprechend 54 dB/ 
Wellenlange. Der Grund fiir dieses Verhalten der 
Dampfungskonstanten ist der gleiche, der zu 
einer fiir alle Frequenzen gleichen Wellenlainge 
an der Resonanzstelle fihrt. 


b) Schnelle der Membran 


Die Schnelle der Membran ist in den Abb. 5 
und 6 eingetragen. Da Schnelle und Amplitude 
iiber die Kreisfrequenz zusammenhangen, kann 
man die Kurven auch als Amplitudenkurven lesen. 
An der Resonanzstelle ist die maximale Ampli- 
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Abb. 7. Phasengeschwindigkeit der Druck- 
welle fiir verschiedene Frequenzen. 


Abb. 8. Dampfungskonstante « fiir 
verschiedene Frequenzen. 
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tude jeweils auf Eins normiert. Man erhalt fir 
die Einhillenden der Amplitudenkurven die be- 
kannten breiten ,,Resonanzkurven“:. Die einge- 
zeichneten punktierten Linien entsprechen wie- 
der errechneten Werten. 

Die Membran ist so bemessen, da fiir eine 
Frequenz von 32 Hz gerade keine Reflexion am 
Ende mehr entsteht. Das Modell umfaf$t einen 
Frequenzbereich von 32:---320 Hz, das sind etwas 
iiber 3 Oktaven (Abb. 13). 

Die Phasenkurven der Membranschnelle bzw: 
-amplitude entsprechen fortschreitenden Wellen. 
Es werden Phasenwinkel von tiber 5z erreicht. 
Die Phasengeschwindigkeit der Welle auf der 
Membran, die kurz als Membranwelle bezeichnet 
werden soll, ist fir die drei Frequenzen 32, 64 
und 160 Hz in Abb. 9 eingezeichnet. Es tritt Dis- 
persion auf, die durch die komplexe Membran- 
impedanz begriindet ist. Bei resultierendem Fede- 
rungswiderstand eilt die Schnelle dem Druck vor- 
aus und bei resultierendem Massenwiderstand 
nach. Durch diese zusidtzliche Phasendrehung 
zwischen Druck und Schnelle ergibt sich fiir die 
Membranwelle eine andere Phasengeschwindig- 
keit als fiir die Druckwelle. Die Abweichung zwi- 
schen beiden ist besonders groB an der Resonanz- 
stelle, da sich dort der Phasenwinkel stark mit 
dem Ort andert. Wie Abb. 9 zeigt, verhalten sich 
die Phasengeschwindigkeiten fiir 160 und 32 Hz 


an der Stelle x=12 cm wie 1 : 4. 
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Abb. 9. Phasengeschwindigkeit der Membranwelle fiir ver- 
schiedene Frequenzen. 


c) Verhdltnis p/v der Membran 


Aus den nach Betrag und Phase gemessenen 
Groen Druck und Schnelle laBt sich das Ver- 
haltnis p/v der Membran bestimmen. Es ent- 
spricht im wesentlichen der Membranimpedanz. 
Die Abb. 10 zeigt den Absolutbetrag und Abb. 11 
den Phasenwinkel fiir die drei Frequenzen. Der 
Betrag fallt mit der Entfernung von der Anre- 
gungsstelle stark ab, woran die Membranbreite 
nicht unwesentlich beteiligt ist. Es war nach 
zl. (3) 

ae? 
of 
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Abb. 10. Betrag der Impedanz p/v fiir verschiedene Fre- 
quenzen. 
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Abb. 11. Phasenwinkel der Impedanz p/v fiir verschiedene 
Frequenzen. 


Da die Flache 6? im Nenner steht, wirkt sie im 
gleichen Sinne wie die mit wachsender Entfer- 
nung kleiner werdende Membranimpedanz Z. Auf 
diese Weise andert sich das Verhialtnis p/v auch 
im Modell um einige GréBenordnungen. Der Ab- 
fall bis zur Resonanzstelle ist so steil, dal} die 
oben beschriebene Verschiebung des Amplituden- 
maximums der Membran von der Resonanzstelle 
zur Anregungsstelle hin nur etwa 0,5---1 cm aus- 
macht. Man kann daher in erster Naherung die 
Membranimpedanz an der Stelle maximaler Aus- 
lenkung als reell ansehen. 


d) Experimentelle Bestimmung der Konstanten der 
Membranimpedanz 


Die Federung wurde statisch durch Auflegen 
von Gewichten und Messung der Membrandurch- 
biegung bestimmt. Sie ist in Abb. 12 fiir die 1 mm 
starke Gummimembran als Funktion der Ent- 
fernung von der Anregungsstelle aufgetragen. 
Ihre GréBenordnung betragt 10-® em/dyn. 

Da der Ort maximaler Membranauslenkung 
als Funktion der Frequenz bekannt ist (Abb. 5, 6 
und 13), kann mittels der Beziehung w-5=1/// MF 
die schwingende Masse berechnet werden. Sie hat 


die GréBenordnung | g und ist in Abb. 12 mit 
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eingezeichnet. Sie ergibt sich proportional zu b*-” 
und entspricht damit ungefahr der Rechnung von 
FLetcHeErR, der ein Anwachsen mit b® angenom- 
men hat. Die Membranmasse, die nur b? propor- 
tional ist, betragt im Mittel nur 1/, der Gesamt- 
masse. 
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Abb. 12. Verlauf des Reibungswiderstandes, der Federung 
und der Gesamtmasse der Membranelemente. 
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Abb. 13. Die Resonanzfrequenz der Membran als Funktion 
des Abstandes von der Anregungsstelle. 


Der Reibungswiderstand ergibt sich aus der 
,.45°-Verstimmung*‘. Aus der Abb. 11 ersieht 
man, da z. B. an der Stelle x= 16 cm eine Fre- 
quenzanderung Af=32 Hz notwendig ist, um 
den Phasenwinkel der Membranimpedanz um 
+ 45° zu veraindern. Man erhalt an dieser Stelle 
(fires = 50 Hz) eine relative Halbwertsbreite Af/f 
= 0,64. Das entsprechende logarithmische Dekre- 
ment betragt 2,1. Es wurde fiir verschiedene 
Stellen der Membran bestimmt, wobei sich Werte 
zwischen 2 und 2,1 ergaben.- Ahnlich den Ver- 
haltnissen im Innenohr ist also das Dekrement 
langs der Membran fast konstant. Unter Zu- 
grundelegung eines Dekrementes von 2 wurde 
der zugehérige Reibungswiderstand mittels der 
bekannten Beziehung fiir die 45°-Verstimmung 


H.-G. DIESTEL: MECHANISCHES MODELL DES INNENOHRES 


429 


berechnet und gleichfalls in Abb. 12 eingezeich- 
net. Der Reibungswiderstand hat die GréBenord- 
nung 10? dyn/em:s*}. 

Mit den Werten nach Abb. 12 wurden die be- 


rechneten Kurven der Abb. 5 und 6 gewonnen. 


e) Eingangsimpedanz des Modells 


Die mechanische Eingangsimpedanz des Mo- 
dells wurde durch Messung der elektrischen Riick- 
wirkung beim Ansetzen des Modells an ein Vibro- 
meter bestimmt [13]. Sie ergab sich in Uberein- 
stimmung mit der Rechnung (Gl. (7)) reell und 
wenig frequenzabhangig p/u = 5,5-103dyn-cm?/ 
em:s, Zwistocki berechnete den Eingangs- 
widerstand des Ohres zu 9-103 dyn: cm~?/em:s"}. 
Reduziert man dieses Ergebnis durch Einsetzen 
des halben Impedanzwertes fiir die Basilarmem- 
bran auf das Einkanalmodell, so erhalt man 
64-103 dyn-cm~/em-s}. Auf den Zusammen- 
hang zwischen diesem und dem am Modell ge- 
messenen Wert wird weiter unten noch einge- 
gangen. 

Bei sehr tiefen Frequenzen ist fir den Ein- 
gangswiderstand des Modells die Federung der 
Anregungsmembran ausschlaggebend, so dab die 
Impedanz dann Federungscharakter annimmt. 


f) Fldachen konstanten Druckes 


Im Gebiet hinter der Anregungsstelle sind die 
Flachen konstanten Druckes zunachst eben und 
stehen auf der Langsachse des Modells senkrecht. 
Bei Anniherung an die Resonanzstelle neigen 
und kriimmen sich die Flachen, wie es in Abb. 14 
z. B. fiir 50 Hz und die Koordinaten x =10 und 
16 cm eingezeichnet ist. Die Krimmung kommt 


Membran { Resonanzstelle 


0 10 16 20 


30 cm 
—— Abstand x 
Abb, 14, Die Flaichen konstanten Druckes (.—+—- ) fiir 


die Frequenz 50 Hz. Oben: schematische Seiten- 
ansicht des Modells; darunter Aufsicht. 


dadurch zustande, da die Membran nicht den 
ganzen Querschnitt des Modells ausfiillt. Es ent- 
steht auch ein Druckgradient in z-Richtung, und 
das Problem wird dreidimensional. Aus diesem 
Grunde und wegen der Querschnittsinderung des 
Kanals kann die zweidimensionale Lisung Gl. (9) 
fiir das Modell nur naherungsweise gelten. Sie 
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gibt fiir den Fall yh =1 eine Phasenverschiebung 
von 50° zwischen dem Druck in Membran- und 
in Bodenn&ahe. 


Experimentell wurden folgende Phasenunter- ° 


schiede festgestellt; z. B. 


f= 50 Hz an der Stelle x=16 cm: 50°, 
f=100 Hz an der Stelle x=10 cm: 30°. 


Die Phasendifferenzen lassen sich auch an 
Hand des Ersatzschaltbildes Abb. 2 verstehen. 
Die Druckmessung in gréBerer Entfernung von 
der Membran, also gewissermaBen auBerhalb der 
mitschwingenden Fliissigkeitsmasse, entspricht 
der Spannungsmessung zwischen den Punkten A 
und B. Geht man mit dem Mikrophon unmittel- 
bar an die Membran heran, so entspricht das 
einer Messung zwischen D und B. Man erhalt 
zwischen beiden Spannungen einen leicht angeb- 
baren Phasenunterschied. Wegen der endlichen 
GréBe des Mikrophons ist es nicht méglich, den 
Druck unmittelbar an der Membran zu messen. 
Man erhalt gute Ubereinstimmung zwischen ge- 
messenen und berechneten Werten mit der An- 
nahme, das das Mikrophon bis in die Mitte der 
mitschwingenden Masse herangefiihrt werden 
kann. Die Messungen entsprechen somit den 
Spannungen zwischen A und B, und zwischen E 
und B. Die Phasenunterschiede zwischen diesen 
bezeichneten Punkten nehmen erst dann gréBere 
Werte an, wenn die Massen- und die Federungs- 
komponente der Membranimpedanz die gleiche 
GréBenordnung haben. Das ist in der Nahe der 
Resonanz der Fall. 


g) Membranauslenkung, Fliissigkeitswirbel 


Die Abb. 15 zeigt vier verschiedene Phasen- 
lagen der Membranauslenkung bei einer Frequenz 
von 32 Hz. Die Membran wurde von der Seite 
aus geringer Héhe photographiert und ist als 


Abb. 15. Vier Phasenlagen der Membranauslenkung; Frequenz 32 Hz. Membran von 
der Seite unter einem Erhebungswinkel von 30° photographiert. Laufrich- 
tung der Welle von links nach rechts. Zwirnsfaden in 5 mm Hoéhe iiber 
der Membran als Bezugslinie. Belichtungszeit 1/1000 s. 
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trapezformige dunkle Flache erkennbar. Sie tragt 
einen weifen Farbstrich langs der Mittellinie, um 
die Wellenbewegung besser sichtbar zu machen. 
Als Bezugslinie wurde ein weiBer Zwirnsfaden in 
5 mm Héhe iiber der Membran ausgespannt. Die 
fortschreitenden Wellen laufen auf den Bildern 
von links nach rechts. Man erkennt das Ampli- 
tudenmaximum und die nach rechts kleiner wer- 
dende Wellenlinge. Im oberen Teil der Bilder 
sind die Schrauben andeutungsweise zu sehen, 
welche die Metalleinsatze, zwischen denen die 
Membran liegt, aufeinanderdriicken. Sie haben 
einen gegenseitigen Abstand von 2 cm. 

G.v. BEKEsy beobachtete im Ohr und in Modell- 
versuchen an der Resonanzstelle einen Fliissig- 
keitswirbel, der in solchem Sinne rotiert, daf die 
Flissigkeit in der Nabe der Membran in Fort- 
schreitungsrichtung der Welle flieBt. ZwisLock1 
beschrieb diesen Effekt durch Beriicksichtigung 
der quadratischen Glieder in den hydrodynami- 
schen Gleichungen. 

Um den Wirbel sichtbar zu machen, wurden 
dem Wasser Aluminiumflitter zugesetzt. Die Be- 
leuchtung der Flitter geschah mit einem Licht- 
biindel geringer Breite, um nur die Teilchen in 
der Mittelebene des Modells aufleuchten zu las- 
sen. Eine Einengung des Fliissigkeitsvolumens 
war nicht notwendig. Die Aufnahmen erfolgten 
senkrecht zum Lichtbiindel. Man erkennt in der 
Abb. 16 die von auBen sichtbaren Schrauben des 
Modells als dunkle Flecke. Die Lange der Bilder 
entspricht 15 cm. An den oberen Bildkanten ist 
der Metalleinsatz als dunkler Streifen zu sehen. 
Die Membran liegt auf den Bildern oberhalb dieser 
Kante. Die Anregungsstelle befindet sich links, 
die Wirbel rotieren im Uhrzeigersinne. Man er- 
kennt, daB der Wirbel fiir héhere Frequenzen 
naiher an der Anregungsstelle liegt. Die Wirbel 
entsprechen Frequenzen von 40 und 80 Hz. 


6. Ahnlichkeit zwischen 
Ohr und Modell 


AbschlieBend sei noch 
kurz die Ahnlichkeit zwi- 
schen den Verhaltnissen 
im Innenohr und dem 
Modell diskutiert. 

Fiir den linearen Ver- 
gréBerungsmaBstab ist das 
Verhaltnis der Membran- 
breiten ausschlaggebend. 
Die Breite der Modell- 
membran Andert sich star- 
ker mit der Entfernung 
als dies bei der Basilar- 
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membran der Fall ist. Der mittlere VergréBe- 
rungsmafstab des Modells betragt 55. 

Beim Vergleich der Stellen = 1,5cem 

(Entfernung vom Steigbiigel im Ohr) 

und x= 16 
(Entfernung von der Anregungsstelle im Modell) 
ergibt sich, da die Werte fiir beide Seiten der zu 
erfiillenden Ahnlichkeitsbedingungen Gl. (11) um 
12°% verschieden sind. Zum Vergleich wurden fiir 
das Ohr die von FLETCHER angegebenen Werte 

M, = 10-531+058, FY, = 109948507 
R, = 0,5 (M,/F,)'" 

herangezogen. Die aus diesen Werten gebildete 
Impedanz der Basilarmembran wurde mit ihrem 
halben Wert eingesetzt,um 
die Verhaltnisse auf das 
Einkanalmodell zu iiber- 
tragen. Der Index 1 soll 
fiir Groen des Ohres, der 
Index 2 fiir Griffen des 
Modells gelten. 

Stimmen die Ahnlich- 
keitsbedingungen in bei- 
den Systemen iiberein, so 
gilt das gleiche fiir den 
Ausdruck yh, worin h eine 
charakteristische Lange 
ist. Der Verlauf der Damp- 
fung und der Phasenge- 
schwindigkeit ist ahnlich, 
das Verhaltnis Wellenlange 
zu Kanaltiefe in beiden 
Fallen gleich. 

Fiir die beiden miteinander verglichenen Stel- 
len des Ohres und des Modells gilt 1, /1, = 55, ferner 
das Frequenzverhiltnis f,/f, = 1/40. Daher stim- 
men die entsprechenden Geschwindigkeiten I, f, = 
0,731, f, der GréBenordnung nach iiberein; ferner 
wegen Gl. (11) auch die Verhaltnisse Druck zu 
Geschwindigkeit. 

So erklart sich die numerische Ubereinstim- 
mung der Phasengeschwindigkeiten in Abb. 7 
mit den von vy. BEKEsy gemessenen Werten. 
G.v.BEKEsy maf in der Nahe des Steigbiigels 
Phasengeschwindigkeiten von 45 m/s, die mit 
wachsender Entfernung bis auf einige m/s ab- 
fielen. 

Ferner findet hier die zunachst zufallig er- 
scheinende Ubereinstimmung der Eingangsimpe- 
danz des Modells mit der von Zwis.ock! berech- 
neten Impedanz ihre Begriindung. 

Im Modell andert sich die Federung der Mem- 
bran nicht so stark mit der Entfernung wie dies 
im Ohr der Fall ist. Daher bleibt die Ahnlichkeit 
zwischen Ohr und Modell nicht iiber die ganze 
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Lange erhalten. Man kénnte daran denken, die 
Federung durch eine allmahlich diinner werdende 
Membran noch starker verainderlich zu machen. 
Will man jedoch das Ohr in allen Teilen ahnlich 
vergréBern, so wiirde eine 55fache VergréBerung 
schon ein 1,8 m langes Modell bedeuten. 

Es wurde daher davon abgesehen, das Modell 
in seiner ganzen Lange dem Ohr 4bhnlich zu 
machen, sondern ein Modell konstruiert, das in 
einem gewissen Bereich die Ahnlichkeitsgesetze 
erfillt und dariiber hinaus qualitativ richtig ist. 


Herrn Prof. Dr. E. Meyer danke ich fiir sein 
forderndes Interesse an der Arbeit und fiir zahl- 


reiche wertvolle Ratschlage. 


(Eingegangen am 1. Oktober 1953.) 


Abb. 16. Aufnahmen des Fliissigkeitswirbels fiir eine Frequenz von 80 Hz (oberes 
Bild) und 40 Hz (unteres Bild). Laufrichtung der Welle von links nach 
rechts. Rotation der Wirbel im Uhrzeigersinn. Membran befindet sich 
einige Millimeter iiber dem oberen Bildrand. Belichtungszeit 25 s. 
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SOME REMARKS ON THE ANALYSIS 
OF MAGNETOSTRICTIVE TRANSDUCERS 


by H. NoprveptT 


Norwegian Defence Research Establishment, Horten (Norway) 


Summary 


The electro-mechanical coupling constant of magnetostrictive vibrating rods is calculated 
from basic physical constants. The formula given is valid for both fundamental and higher 
harmonics. The calculated values are compared with directly measured coupling constants. Next, 
the electric analogue of rods is given, and it is shown how a graphical analysis of admittance and 
impedance diagrams gives convenient rules for evaluating the relevant quality factors. Actual 
measurements on a particular specimen are included. 


Sommaire 


On calcule, 4 partir des constantes physiques fondamentales, le coefficient de couplage électro- 
magnétique des tiges vibrant par magnétostriction. La formule établie est valable 4 la fois pour 
la fréquence fondamentale et pour les harmoniques d’ordre élevé. On compare les coefficients de 
couplage ainsi calculés avec ceux mesurés directement. On indique ensuite l’équivalent électrique 
d’une tige et on montre qu'une étude graphique des diagrammes d’admittance et d’impédance 
permet d’établir des régles de détermination des facteurs affectant la qualité. On donne enfin les 
résultats de mesures faites sur une éprouvette. 


Zusammenfassung 


Der elektro-mechanische Kopplungsfaktor fiir stangenfoérmige Magnetostriktionsschwinger 
wird aus den fundamentalen physikalischen Konstanten berechnet. Die angegebene Formel gilt 
sowohl fiir die Grundschwingung als auch fiir hdhere Harmonische. Die berechneten Werte werden 
mit gemessenen verglichen. Weiterhin wird ein elektrisches Ersatzschaltbild fiir die Schwinger 
angegeben, und es wird gezeigt, wie man durch graphische Auswertung der Admittanz- und 
Impedanzkurven die Giite bestimmen kann. Fiir ein spezielles Beispiel sind Messungen angegeben. 


1. Introduction and phase relations along their lengths, and they 
have an external magnetic resistance which can- 
not be neglected. We shall see that more general 
expressions for the electro-mechanical coupling 
may be given for these cases. 

The electric network analogues which can be 
substituted for the magnetostrictive vibrating 
system are of great value in design and measure- 
ment of transducers. It is felt, however, that full 
use of these analogues has not been made, and 
that principles and methods used in other fields 
of electrical engineering may be adopted to in- 
crease their usefulness. In the following we will 
show that graphical methods, applied to the 
admittance and impedance loci, will yield the 
determining factors of the transducer, more 
quickly and easily than an analysis of the 
algebraic equations. 


Although the magnetostrictive transducer has 
achieved a considerable importance in practice, 
particularly in underwater acoustics, it cannot 
be said that the theory is developed to the same 
extent as is the case with the piezoelectric trans- 
ducer. In most treatments en magnetostrictive 
transducers, the paper by BurreRWORTH and 
SmiTH [1] is usually referred to when dealing 
with the underlying principles. This paper con- 
sidered the ring-shaped transducer, excited in its 
fundamental radial mode, and an equivalent 
electric analogue was worked out, together with 
an expression for the electromechanic coupling 
coefficient, from the differential equations of the 
system. The ring-shaped transducer was charac- 
terised by having a closed magnetic circuit, and 
by having all elements round a periphery moving 
with the same amplitude and in the same phase. 
The electro-mechanical coupling coefficient was, 
in this particular case, expressed by the physical The electro-mechanical coupling coefficient 
constants of the material only. Since the appear- expresses to what extent the vibrating mechano- 
ance of the paper mentioned, several transducer acoustic system is coupled by the magnetostric- 
types have been developed, which show neither — tive effect to the attached electric network, and 
of the two characteristics. Longitudinally vibrat- is the most important quality factor of the 
ing rods, for instance, exhibit varying amplitude transducer. The more the coupling factor exceeds 


2. Electro-mechanical coupling relations 
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its minimum value, determined by the losses in 
the electric and mechano-acoustic parts, the 
more efficient is the energy conversion of the 
transducer. The bandwidth over which it may 
be used efficiently, depends on the ratio of the 
actual and the minimum coupling value. An 
evaluation of this quality factor is therefore of 
interest. 

This evaluation is usually carried out by 
determining the electric analogue and then giving 
the coupling coefficient as the square root of the 
ratio of the two inductances appearing in the net- 
work. The coefficient can, however, be found di- 
rectly, if we consider the energy conditions in the 
combined electro-mechanical system. The electro- 
mechanic coupling coefficient may be defined as 
the square root of the ratio of the magnetic ener- 
gy stored in the system, to the stored mechanical 
energy present. Applying this point of view first 
to the ringshaped transducer, vibrating in its fun- 
damental mode, we obtain the stored mechanical 
energy (when it is all in potential form) in a 
circumferential element of length ds as: 

rn. = ABA (=)'as 
2 s 
where E is the elastic modulus of the material, 
A the cross-section of the ring, € the vibrational 
amplitude (with reference to change of circum- 
ference) and s the circumference corresponding 
to the average diameter of the ring. 

The magnetically stored energy in the same 

element is given by: 


where « is the inverse magnetostrictive constant 
connecting the developed magnetic field H with 
the relative change of length, according to the 
equation: 

E 

H=«—., 

s 

while « is the permeability of the material. 
Evaluating the coupling coefficient from the 


energy ratio, we obtain: 


sees =) eu  |/Kap 

ahi 4nE om 
where A= «/4z, the magnetostrictive constant. 
The expression for the coupling coefficient is 
identical with that evaluated by BurreERworTH 
and Smira. 

Dealing with vibrating rods, the conditions 
are somewhat different. We will consider a 
magnetostrictive rod of length / in a medium, 
excited in a longitudinal vibration by, for ins- 


k= 
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tance, an incoming sound wave of angular fre- 
quency w. The displacements of the elements 
along the rod can then be written as: 


— = cos = (l—x) 


where c is the sound velocity in the rod, and x is 
the distance from the excited end, the other end 
being free. We will for the sake of simplicity 
assume that the rod has a cross-sectional area of 
1 cm?, and we will further assume that there is a 
steady polarising field By along the length of the 
rod. The effect of this field is to make the connec- 
tion between the superimposed alternating field 
and the mechanical displacements unambigous. 
It does not enter into the subsequent evaluation. 

Since the alternating field developed by the 
vibration at any point along the length of the rod 
is given by d(4H) =x (dé/dx) dx, the resulting 
field is determined by the relative displacement 
of the two ends. The difference in displacements 
at x=0 and x= 1 is easily found to be: 


6) =F (1 cos = i) ; 
Consequently, the developed induction is 
given by: 


A Benaee: (1—cos 1) 
l c 


where « now is the effective permeability of the 
rod and not that of the material. 

The expression for the alternating magnetic 
field has maxima for m=(2n+1)zc/l, i.e. for 
the fundamental and all odd harmonic modes. 
For even harmonic modes, no magnetic field will 
be developed by the vibrations. 

Calculating the magnetically stored energy in 
the rod in the odd modes, we obtain: 


tt pn €2 
| Aare ric HEA ee Say hes 


while the mechanical energy stored in potential 
form is found from: 


dé \? 1 ])2 72 
x 4, 


Again evaluating the coupling factor from 
these energy ratios, we arrive at the formula: 
2/2 / KA 
(Qn+1)x¥ E - 
This formula contains the expression for the 
coupling factor of a vibrating ring, but has an 
additional factor due to the variation of displace- 
ment along the direction of vibration. The 
reducing effect of this factor will, in the funda- 
mental mode, be (2|/2)/7= 0.9, but will in the 
higher modes have a much more serious effect. 


Wek = 
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How the difference in coupling values at the 
even and odd modes may be used to separate the 
A—E effect appearing in some ferromagnetic 
materials, from their magnetostrictive effect, has 
been shown in an earlier paper [2]. 

The evaluated formula will now be verified by 
experiments. A rod of length 10 cm and with a 
cross-section 11cm? is built up of annealed 
nickel stampings by using some sort of cement. 
The physical constants contained in the coupling 
formula, i. e. the permeability . and the magneto- 
strictive constant A, are then measured, and the 
value of the coupling constant found from the 
formula. Afterwards, the coupling factor will be 
measured directly by examining the impedance 
conditions at the terminals of a surrounding coil. 

Since we have assumed a polarising field to be 
present, « and / will be incremental values. 
Although incremental values for the permea- 
bility of nickel are given in various handbooks, 
we prefer to measure this constant on the spe- 
cimen in question, since magnetic constants have 
a tendency to vary considerably with heat treat- 
ment and cold working. 

Confining our attention to the coupling condi- 
tions at the fundamental frequency of the rod, 
the incremental permeability should be measured 
at this particular frequency. The vibrations of 
the resonant rod will, however, influence the 
magnetic measurements, unless a complete 
mechanical clamping of the rod is obtained. This 
being difficult, it is more practical to do the 
measurement for two frequencies—one a little 
below and one above. The values for the incre- 
mental permeability of the rod given in Fig. 1 
are obtained by taking the average of two such 
measurements. 


4 ~ 10° qouss 6 


Fig. 1. Incremental permeability of nickel rod versus 


polarising induction. 


Curves for the variation of length with 
magnetic field or induction are given in several 
books on magnetism, but values for the incre- 
mental change of length with varying polarising 
field are scarce and none seem to be available for 
nickel. The assumption that the incremental 
value of A is equal to the slope of the static curve 
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does not seem to be justified. It appears that a 
difference exists and the constant has to be 
measured directly. This is done by using a C. E. 


’ Johansons Mikrokator which is able to measure 


mechanical displacements less than 10 mm. By 
superposing a small incremental induction on 
polarising fields of varying magnitude, the curve 
in Fig. 2 is obtained. 


12 


~10 °gauss 


8 


A(Al/L) 
4B 


0) 2 4x 10° gauss 6 


Fig. 2. Relative incremental deformation of nickel rod, 
versus polarising induction (A= E-A (Al/l)/AB). 


As is shown, there appears to be a saturation 
value for the incremental magnetostrictive con- 
stant, something that would be rather difficult to 
discern from the static curve. Usually the static 
change of length is supposed to vary quadrati- 
cally with the magnetic field, giving a linear 
increase in the incremental value with increasing 
polarising fields. 


20 


1, 


0 2 410° gauss 6 


Fig. 3. Coupling coefficient of nickel rod versus polarising 
induction (© are directly measured values). 


The third factor in the formula for the coupling 
coefficient, the elasticity modulus, may safely be 
taken out of handbooks as being 2.1 10" 
dynes/cm?”, and inserting the three constants, we 
obtain the curve in Fig. 3 for the coupling coeffi- 
cient of the rod as a function of the polarising 
field. This curve will now be checked by direct 
measurements at the terminals of a coil surround- 
ing the rod. We will then use the method given 
by Rounp [3]. In short, this is based on the fact 
that the impedance of a magnetostrictive trans- 
ducer, tuned with a capacitance at its resonant 
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frequency, exhibits three resistive positions 
around this frequency. One of these positions will 
be at the actual resonance frequency f), the other 
two will be displaced by +Af respectively. The 
coupling coefficient is then quite simply given by 
2Af/f,. Carrying out this measurement on the 
rod, we obtain the values indicated in Fig. 3. It 
appears that good agreement exists between the 
two methods of evaluating the coupling coeffi- 
cient. 


3. Electric analogues 


Having worked out the general expression for 
the coupling constant, it is easy to determine the 
electric network corresponding to the magneto- 
strictive rod. It is apparent from the preceding 
section that excitation of the rod in higher 
harmonics is not to be recommended if good 
conversion efficiency is required. We will there- 
fore confine ourselves to the fundamental mode, 
although analogues for the higher modes may be 
evaluated in the same way. 

Still considering longitudinal vibrations only, 
the impedance offered by a rod to an incoming 
sound wave, is given by: 
velocity . 1 


a) 
= cot—il 
pressure J oc c 


Vices = 


(o being the density of the rod) using the so 
called mobility or electro-magnetic system of 
equivalents. This system must be used when 
dealing with magnetostrictively coupled trans- 
ducers, as pointed out by Mixes [4]. 

Considering first the non-magnetostrictive rod, 
it is easy to show that Z,,.., around the funda- 
mental resonance is closely approximated by the 
impedance of an electric parallel circuit, where 
the capacitance is given by Cyn =—ol/2, i.e. 
half the weight of the rod, and the inductance by 
Lnech = 21/x? E, i. e. 2/x* times the static compli- 
ance of the rod. In this network, the vibrational 
energy will alternately be stored in elastic and 
motional form. 

Assuming now the rod to be magnetostrictive, 
there must be a storage element in the analogue 
for the magnetic energy k? Wiech. We obtain a 
network expressing this physical condition by 
adding an inductance k? Lj, to the purely 
“mechanical” inductance evaluated.Consequently, 
we get the network given in Fig. 4a, as represent- 


(a) (b) 


Fig. 4. Electric analogues of nickel rod. 
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ing the vibrating magnetostrictive rod. It will be 
noticed that no mention has been made of any 
surrounding coil in this evaluation of an electric 
analogue. The rodis, initself, an electro-mechanical 
system where both mechanical and _ electro- 
magnetic quantities must appear, if an analogue 
is worked out. 

The rod has up to now been assumed loss free. 
In practical cases, losses occur due to inner 
friction in the material and to radiation of sound 
into the surrounding medium. These losses may 
be represented by a resistance across the ter- 
minals of the condenser in the analogue. The 
value of this resistance will, in addition to the 
inner friction in the material, depend on the 
acoustic impedance of the medium, and on the 
area of the radiating end in relation to the wave- 
length of the transmitted sound. 

If no magnetic leakage is assumed, a coil 
surrounding the rod will appear in the analogue 
as an inductance closely coupled to the magneto- 
strictive inductance k? Lynn. The value of this 
purely electric inductance is worked out from 
the standard formula for inductances, i.e. 
Leai=42N*u/l (IN being the number of turns), 
and the losses due to eddy currents may be 
taken into account by adding a parallel resist- 
ance across the terminals. The complete network 
equivalent for a magnetostrictive rod transducer 
finally takes the form given in Fig. 4b, where all 
components are given by basic physical constants. 


4, Analysis of admittance and 
impedance diagrams 


The analogue given in Fig. 4b is a little incon- 
venient to analyse graphically, due to the parallel 
combination of the capacitance and resistance. 
Since most magnetostrictive transducers are 
highly resonant, we may without any great error 
replace the parallel combination by a series RC 
branch. If we, in addition, replace the electro- 
mechanical transformer by its equivalent, we 
obtain a network of the type given in Fig. 4c, 
which is more suitable for practical purposes. 

Using well-known methods, as outlined for 
instance by M.1.T. [5], the admittance diagram 
of Fig. 4c is evaluated. Inverting the line repre- 
senting the impedance of the R’ L’ C’ combination 
[Z’ = R'+j (wL'—1/C’)], we obtain the circle 
around M in Fig. 5, where o is the origin. The 
admittance of the electric branch (R/,L;,) may 
be added by moving the origin to O, where OA = 
l/oLi; and Ao=1/Ri,. This latter admittance 
will vary with frequency, although this variation 
is small compared to that taking place in the 
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mechanical branch, if the rod is vibrating in air. 
If the rod is in a medium of less acoustic mobility, 
however, this variation must be taken into 


account. As is shown for instance by BERANEK: 


[6, p. 353], the contribution of the mechanical 
branch may be separated, by assuming the 
electric admittance to vary regularly through the 
resonance point. The motional impedance circle 
can then be constructed by plotting the vectorial 
difference between the admittance actually 
measured, and the value of the electric admit- 
tance at each particular frequency. When this 
operation is carried out, we obtain an admittance 
diagram of the type given in Fig. 5, such that we 
may confine ourselves to a treatment of this. 


0 ee QU 


Fig. 5. Evaluation of admittance diagram for magneto- 
strictive transducer. 


The point b in Fig. 5 obviously corresponds to 
the resonance frequency of the rod, since the 
mechanical capacitance and inductance cancel at 
this point. At the point a and ¢ the reactive 
component of the motional admittance is equal 
to the resistive, and we obtain the expression for 
the mechanical Q-value of the rod as: 


(alee 
seh PS fac 

The Q-value of the coil itself is of course given 
by the ratio of the electric resistance to the 
inductance. 

The electric energy absorbed by the network 
in Fig. 5 for any frequency around the resonance 
point is proportional to the distance from the 
ordinate axis to the point on the circle at the 
particular frequency. The ordinate axis may 
therefore be considered as the locus of the points 
in the plane with zero absorbed energy. 

Assuming that the vector Oo is the value of 
the electric admittance at resonance, and that we 
consider a limited band around this frequency, 
the energy absorbed by the mechanical branch 
in the network is proportional to the distance 
from the line m—m to points on the admittance 
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circle for the various frequencies. The line m—m 
can consequently be looked upon as the locus for 
the points in the plane with zero mechanically 
converted energy. The idea of loci for various 
energy conditions will be useful when we later 
consider the impedance diagram. 

The maximum electro-mechanical conversion 
efficiency occurs at the resonance point b, and 
can be expressed as the ratio: 

Hreck OB Ab: 

The admittance diagram is convenient to 
analyse. When, in spite of this, the impedance 
diagram is usually adopted in practice, it may 
be because the impedance concept is more 
familiar to engineers, and because the majority 
of bridges are of the impedance type. 

The principle of analysis can, however, also be 
applied to the impedance locus and, although the 
evaluation is a little more complicated, the 
quality parameters of the transducer are just as 
easily taken out of the finished diagram. 

The impedance of the circuit in Fig. 4¢ can be 
written as: Zz 


2 = Li 
a nck + La 


ofa 
Introducing Z/;=|Z21| e (3 fou) where 9 = 
tan! (w L;1/R2i) is the loss angle of the coil, we 
can write: 
Zire” 
Linoch Liat” 


The locus of the vector expressed by the 
second term is given by the circle around M, 
inclined at an angle 2 9 to the horizontal in Fig. 6. 
The impedance of the electric branch is again 
added by moving the origin from o to O, where 


L= Leia 


Oo is the impedance of Ri; and L/,. The same 
considerations must be applied regarding the 
variation in Z,, as was necessary when dealing 
with the admittance diagram. 


> 
rs) 
= 
& 
3 
= 
lu 


Fig. 6. Evaluation of impedance diagram for magneto- 
strictive transducer, 
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We will now show how the relevant points in 
the diagram may be found. Assuming that the 
two diagrams in Figs 5 and 6 represent the same 
circuit, one diagram must be the inversion of the 
other, by definition, i. e. the impedance diagram 
can be inverted into the admittance diagram, 
and vice versa. By choosing the constant of inver- 
sion properly, the two diagrams may be made to 
coincide. Consequently, the point o in Fig. 5 
corresponds to point o in Fig. 6, the point a in 
Fig.5 to the point a in Fig.6, etc. We can there- 
fore find the point corresponding to resonance 
in the impedance diagram by drawing the line 
Ob’, intersecting the circle at b. The inverted 
quadrantal positions can be found in the same 
way by drawing the lines Oa’ and Oc’. The 
mechanical Q-value is consequently given by the 
same expression as in the admittance case. 

The point on the circle corresponding to the 
resonance of the mechanical branch in series with 
the electric branch is more easily found from the 
impedance diagram. At this frequency the reac- 
tive part of Zirecn + Zi; is zero and the second 
term in the impedance expression will have 
maximum magnitude, i.e. the corresponding 
point on the circle will be d. This point can then 
be found in the admittance diagram by drawing 
the line Od’ in Fig. 5. 

Since the coupling constant can be expressed 
by the ratio of the resonance frequency of the 
mechanical branch to the resonance frequency of 
the mechanical and electric branches in series, 
we have found a method of evaluating this 
constant. 


We obtain: ,£7~ V 1—(fa/fs)? 
valid for both the admittance and the impedance 
diagram. 

As before, we can find loci for the various 
energy conditions. The ordinate is obviously the 
locus for the points of zero supplied electric 
energy. Similarly, a line representing the points 
in the plane where no energy is absorbed by the 
electric branch can be evaluated. This evaluation 
is carried out by Bracstap and SKANCKE [7, p. 
46] and we will confine ourselves to giving the 
method of constructing the line. 

The two tangents to the circle from the origin 
O are constructed, and the line m—m is drawn 
through the points of contact. The line OT is 
halved, and a perpendicular n=n is erected at the 
midpoint. This line represents the points in the 
plane where no energy is dissipated in the coil. 

The line for zero absorbed energy in the mech- 
anical branch must obviously go through the inter- 
section point S between the supplied energy line 
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and the electric loss line. It must also go through 
the point 0, since at that point no energy is 
supplied to the mechanic branch. 

Still referring to the work by Bracstap and 
SKANCKE, it can be proved that a linear scale for 
the electro-mechanical conversion efficiency may 
be constructed. Drawing the line h—h intersecting 
the lines n—n and g-g, we can find the efficiency 
at any point on the circle by drawing a line from 
that point through S until intersection with h—-h. 
The point of 100% efficiency is the intersection 
point between h—h and n-n, while the zero point 
on the scale is the intersection point between h—h 
and g-g. The maximum efficiency occurs at the 
resonance point b on the circle, something which 
is expected and which can be proved analytically. 

The efficiencies read on the scale give the pro- 
portion of the supplied electric energy which is 
converted into vibrational energy. How this 
energy is divided between the internal losses 
occuring in the rod and useful radiated energy 
may be found by considering the connection bet- 
ween the resistance R’ in Fig. 4c and the dia- 
meter of the admittance or impedance circle. If 
the rod is vibrating in air, the radiation losses are 
negligible and the resistance, which may be 
denoted by Ri, is almost entirely due to inner 
losses. Submerging the rod in water or another 
medium with small acoustic mobility will in- 
crease the resistance in the network by an amount 
Ri,. The mechano-acoustic efficiency can then be 
written as: Ri 

ee RR 

In the admittance diagram, the diameter of 
the circle in air is given by D,=1/Ri, and of 
that in the water by D.=1/(Ri+ R..), so that 
the efficiency is quite simply given by: 

Da— Dw 
Tima = —y + 

In the impedance diagram Da =) Z/,|?/(Ri+ Ro) 
while Dw = |Z¢i|?/(Ra + Ri. + Roi), where Ré is the 
equivalent series resistance of R/). 

Evaluating, we obtain: 

1 — (Dw/ Daz) 
Tah Z 


Since the second term in the denominator is 


Uf] m-a — 


very small in practical cases, the expression for 
the efficiency takes the same form as in the admit- 
tance case. 


5. Example 


Using the rod described in Section 2, we have 
measured the impedance diagram when the 
surrounding coil has 60 turns and the polarising 
field is adjusted to 3350 gauss. A Birmingham 
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Sound Reproducer frequency generator with 
incremental frequency scale has been used 
together with a General Radio Impedance 


Bridge, Type 916-AL and the diagram in Fig. 7 © 


obtained. There are two circles in the diagram, 
representing the fundamental and third harmonic 
resonance. 


Fig. 7. Measured impedance diagram of nickel rod. 


Utilising the methods evaluated earlier, we 
first locate the fundamental resonance frequency 
fo as being 24350 c/s, and the quadrantal posi- 
tions f. and f- as 24450 ¢/s and 24270 c/s respec- 
tively. The mechanical Q-value of the rod vibrat- 
ing in air is consequently Qmech = fo/(fa—fc) =135. 

The frequency corresponding to the resonance 
of the mechanical and electric branches in series 
is found to be fa= 24140 c/s, and the coupling 
coefficient is consequently given by: 

k =V1—(falfs)? x 100% = 13.2% 
as compared to 15% taken out of the earlier 
given curve in Fig. 3. 

An efficiency scale is constructed and we read 
a maximum electro-mechanic conversion effi- 
ciency of approx. 95% at the resonance point. 

From the measurements we can determine the 
components in the electric analogue. We have: 


Ro + Ra =| Zel?/D. = 287/58 135 ob 
from which: 
Ri=13.5—4.7 = 8.8 ohm. 


The inductance in the mechanical branch is given 
by: 
— Qmech Ra 135 x 8.8 x 10° 


san Sha A) EU igual 
from which the capacitance is found as: 
Geel ut! Bupa Jl oa 
Si ree By (6.28 x 24350)? 7.75 x 10-8 
= 0.0055 uF. 
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Expressing the electric inductance as kL’, we 
obtain when we use the coupling value from the 
curve in Fig.3 (which is probably the more 
accurate) : 


Liy= 0.0225 x 7.15 x 10° = 0. Lisa 


a value which agrees quite well with that taken 
directly out of the impedance diagram. 

In this particular case, when the vibrating 
body is of a simple geometrical shape, it should 
be possible to determine the components in the 
network directly from the basic formulas given 
earlier. 


We have: 
ee 2 Jl tage 10 
mech 2 AE 7? 1x21 102 
= 0.965 x 10-2 CES 
and 


Cmnech = 5 Or AL = > x 8.9 x 1x 10 = 44.5 CGS. 


The electro-mechanical transformation ratio is 
found from Fig. 4b as: 
mNAE  260x1x2.1x10¥ 
41) 4x 10x 1.35 x 104 
= 25% 10", 
The elements appearing in the electric ana- 
logue should then be: 
DE! = T? Lincs = (23.10%)? 0.905 3a ree 
= 5.1.x 108.GCS =:S2 i makes 
C’ Crech 44.5 
frep is (23 x 108)? 
= 8.7 x 10-7 = 0.0087 
The discrepances between these values and 
those evaluated from the impedance diagram 
may be due to errors in the determination of the 
magnetostrictive constant 4 and the mechanical 


Q-value. Also, the stray inductance has not been 
taken into account. 


T= 
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NOISE LEVELS UNDERNEATH SOME CIVILIAN AIRCRAFT 
SHORTLY AFTER TAKE-OFF 


by P. H. Parkin and H. J. Purxis 


Building Research Station, Garston, Herts, England 


Summary 


The maximum sound pressure levels (in octave bands) have been measured underneath several 
types of civilian aircraft shortly after take-off. The height of the aircraft varied from 15:::90 m; 
the results have been corrected to a standard height of 38 m. 


Sommaire 


On a mesuré, dans les bandes de fréquences d’une octave, les niveaux maxima de la pression 
sonore sous différents types d’avions civils, peu aprés leur décollage. L’altitude des avions variait 
de 15---90 m; on a ramené les résultats 4 une altitude uniforme de 38 m. 


Zusammenfassung 


Der maximale Schalldruckpegel (in Oktavbindern) unterhalb verschiedener Typen von Zivil- 
luftfahrzeugen kurz nach dem Start wurde gemessen. Wegen der verschiedenen Flughéhe (15-:-90 m) 
wurden die Ergebnisse unter Annahme eines 1/r?-Gesetzes auf eine Bezugshéhe von 38 m bezogen. 


1. Introduction 


As part of an investigation into the sound insu- 
lation of buildings at London Airport, measure- 
ments were made of the noise from civilian air- 
craft. Some of these measurements were of the 
sound pressure levels in octave bands under the 
aircraft shortly after take-off; as such results do 
not appear to have been published before (except 
for a few individual aeroplanes), and as they may 
have some general interest, a brief summary of 
them is given in this paper although, as will be 
obvious, they by no means constitute a complete 
investigation. 


2. Technique 


A Standard Telephones and Cables moving- 
coil microphone type 4021, protected by a double 
wind-shield, was used at a height of 1.5 m (5 feet) 
above the ground just off the end and on the 
centre-line of the main runway. As the aircraft 
took off they flew over the microphone at heights 
varying from 15 m (50 feet) to about 90 m (300 
feet); their height when over the microphone was 
measured to the nearest 7.5 m (25 feet) using a 
simple theodolite. The noise was recorded on 
magnetic tapeand subsequently analysed through 
octave filters (nominal pass-bands, 37:---75 ¢/s, 
75-150 ¢/s, ete.) into a logarithmic level recorder, 
corrected to read approximately mean rectified 
levels. The results given are derived from the 
maximum values reached by smooth lines drawn 
through the traces of the level recorder. 


It was desirable to correct the results to a 
standard height of the aircraft above ground and 
although there were not sufficient measurements 
to be able to be precise about the accuracy, the 
following procedure appears reasonable. Of the 
four types of aircraft of which at least ten meas- 
urements were made, two types (the “Eliza- 
bethan” and the “Viscount”) showed a mean 
variance (derived from the variances at parti- 
cular heights for each type) of the order of 2 to 
3 dB? in most of the octave bands, while the two 
other types (the “Convair” and the “DC-4’’) 
showed a mean variance three of four times as 
great. (It may be only a coincidence that the 
“Elizabethan” and the “Viscount” are compara- 
tively new aircraft—first made in 1951—while 
the “Convair” and the “DC-4” are comparatively 
old—first made about 1947.) The variation of 
sound pressure level with height above ground of 
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Fig. 1. Variation with height of sound pressure levels 
(300---600 c/s band), Elizabethan e and Viscount 
Aircraft x. 
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Table of results 


Engines 


Maximum loaded 
weight in kg (Ib.) 


Total maxi- 
mum 
horse-power 
at take-off 


Type 


meas- 
ured 


Number 


PISTON ENGINES 


DC-3 2 2 400 11441 (25 200) 5 
Viking 2 3 380 15 422 (34000) 1 
Convair 2 4 800 18972 (41790)| 15 
Elizabethan 2 5 200 23.600 (52000)} 11 
(Ambassador) 

DC-4 4 5 800 33142 (73000)| 13 
York 4 6 480 29500 (65 000) 1 
Argonaut 4 7040 37300 (82 300) 4 
(Canadair) 

Hermes 4 8 400 39 100 (86000)} 3 
DC-6 4 9 600 44.130 (97 200) 1 
DC-6B 4 9 600 48 125 (106 000) 2 
Constellation 4 10 000 48 580 (107 000) q 
Strato-Cruiser | 4 14.000 76 195 (145 800) 2 


PROPJET ENGINES 


4 | 5600 H.P. 
plus 660 kg 
(1 460 Ib.) 


static thrust 


TURBO-JET ENGINES 
| 4 | 9100 kg 


Viscount 23 814 (52500); 10 


Comet 45540 (105000)| 1 
(20 000 Ib.) 


static thrust 


the “Elizabethan” and the “Viscount” agreed 
reasonably well in most of the octave bands with 
the inverse square law (e.g. Fig. 1) and it was 
therefore thought justified to correct all the re- 
sults to a standard height, 38 m (125 feet). A 
further indication of the validity of this proce- 
dure was that the variances, now calculated as- 
suming the inverse square law, were not generally 


changed. 
3. Results 


The results are given in the Table, which shows 
the maximum sound pressure level reached in 
each octave band during passage of the aircraft 
overhead shortly after take-off. All data are cor- 
rected to an aircraft height, at the overhead posi- 
tion, of 38m (125 feet) assuming the inverse 
square law. Where measurements were made on 
more than one aircraft of the same type the fig- 


Means of maximum sound pressure levels in dB above 0.0002 dynes/cm?; 
(standard deviation); 38 m (125 ft.) height 


600--- | 1200-+- | 2400--+ | 4800- - - 
1200 c/s | 2400 c/s | 4800 ¢/s | 9600 c/s 


99 | 106 | 104 | 101 | 100 | 98 | 93 | 80 
101 | 103 | 102 | 103 | 104 | 103 | 100 | 94 

98 | 108 | 104 | 102 | 101 | 99 | 93 | 82 
(1.5) | (4.0) | (2.5) | (3.0) | (3.0) | (2.5) | (3.5) | (4.5) 

95 | 104 | 99 | 94 | 92 | 86 | Shemerge 
(2.0) | (2.0) | (2.0) | (1.0) | (1.0) | (1.5) | (1.5) | (1.5) 
103 | 112 | 107 | 104 | 106 | 104 | 102 | 95 
(2.5) | (3.5) | (3.5) |-(40) | (2.5) | (2.5) | (2.5) | (3.0) 

98 | 108 | 110 | 108°] 108 | 109 | 109 | 104 
100 | 116 | 109 | 102 | 105 | 107 | 106 | 101 
103 | 112 | 106 | 106 | 109 | 104 | 103 | 99 
104 | 111 | 112 | 109 | 108 | 108 | 105 | 100 
104 | 116 | 113 | 110 | 110 | 109 | 108 | 102 
105 | 113 | 110 | 110 | 112 | 109 | 108 | 99 
111 | 117 | 114 | 108 | 107 | 108 | 106 | 97 
101 | 105 | 99 | 96 | 103 | 102 | 101 | 100 
(2.0) | (3.5) | (1.5) | (1.0) | (1.5) | (1.5) | (2.5) | (2.0) 
log | 119 | 124 | 123 | 122 | 119 | 114 | 108 


ures given are the means of these corrected values. 
The standard deviations from these means, where 
calculated, are given to the nearest 0.5 dB in 
brackets. It should be mentioned that the maxima 
in the various octave bands do not necessarily 
occur at the same instant of time, due to the 
directionality of the source. The details of the 
aircraft are taken from ‘“‘Jane’s All the World’s 
Aircraft”, (Sampson Low, London). 
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Sommaire 


Exposé des différentes maniéres de considérer la pression de radiation et de la calculer. Gran- 
deur scalaire ou grandeur vectorielle? L’auteur introduit le vecteur radiant égal au produit de 
la surpression par la vitesse de déplacement d’une tranche d’onde plane et calcule la pression 
de radiation sur un obstacle en considérant séparément les ondes incidentes et les ondes réfléchies. 
Il étend son calcul au cas des ondes planes progressant dans un milieu absorbant et A celui des 
ondes sphériques. 


Zusammenfassung 


Es werden die verschiedenen Moglichkeiten, den Schallstrahlungsdruck zu betrachten, zu- 
sammengestellt. Ist er ein Skalar oder ein Vektor ? Der Autor fiihrt einen ,,Strahlungsvektor“ ein, 
der gleich dem Produkt aus dem Uberdruck und der Verschiebungsgeschwindigkeit einer ,,Scheibe“ 
aus einer ebenen Welle ist, und berechnet den Strahlungsdruck an einem Hindernis, indem er 
die einfallenden und refiektierten Wellen getrennt beriicksichtigt. Die Berechnung wird auch auf 
den Fall einer ebenen Welle in einem absorbierenden Medium und auf Kugelwellen ausgedehnt. 


Summary 


An account of the different ways of considering radiation pressure and of calculating it. Is it a 
scalar or vectorial quantity ? The author introduces a ‘radiation vector’ equal to the product of 
excess pressure and velocity of displacement of a section of a plane wave and calculates the 
pressure of radiation on an obstacle, considering incident and reflected waves separately. He 
extends his calculation to the case of progressive plane waves in an absorbent medium and to 


that of spherical waves. 
Introduction 


La notion de pression de radiation qui a été 
introduite par MAXxwELt il y a trois quarts de 
siécle et qui s’est naturellement étendue a toutes 
les ondes de nature électro-magnétique s’étend 
également au domaine des ondes acoustiques. 
Cependant dans beaucoup de traités, on ne trouve 
rien a son sujet, et, dans les autres, cette pression 
est congue différemment suivant les auteurs: en 
particulier, c’est sous langle vectoriel qu’on la 
trouve envisagée dans les travaux de L. BriL- 
LOUIN [1], tandis qu’elle est essentiellement con- 
sidérée comme une grandeur scalaire dans les 
études de Y. Rocarp [2]. D’autre part, les mé- 
thodes de calcul sont trés variées et généralement 
on se contente de déterminer une valeur moyenne. 


1. Critique d’un mode classique de calcul 


Avant d’exposer notre propre conception et de 
proposer une solution particuliérement simple 
du probléme, nous voudrions reproduire l'un des 
modes de calcul les plus répandu afin d’en faire 
la critique et d’en discuter la validité. 


a) Soient done des ondes planes progressives 
se propageant dans la direction des x croissants 


et définies par le déplacement élémentaire X 
d’une tranche d’abscisse x au repos 


X= X)sin(wt—qx). 
La vitesse de cette tranche a pour expression: 
u= 0X /0t = X) w cos(wt—qx) 
et la surpression Ap correspondante: 
Ap =—P @X |éx = PX,qcos(wt—qzx), 


P étant Vinverse du coefficient de compressibilité. 
Si la propagation se fait dans un gaz d’une fagon 
adiabatique, P a pour expression yp, p étant la 
pression du gaz au repos et y étant C,/C,. La sur- 
pression en un point d’abscisse x ne serait pas 
Ap, celle que nous venons de définir, car l’abscisse 
de la tranche considérée n’est plus x, mais (x + X). 
Elle aurait pour expression: 


Ap, = PXyqcos[wt—q(x+ X)] 
soit en développant: 
Ap, = P Xyq[cos(wt—qx) + qX sin(wt—qx) + -] 
ou encore: 
Ap, = PX oq cos(wt—qx) + 
+ PX? q*sin?(wt—qx)+ --]. 


4.42 


Le premier terme correspond a la surpression 
«normale». Le second serait relatif a la pression 
de radiation: 


Apr= PX} q@ sin? (wt—qx) 


avec pour valeur moyenne P X@ q?/2. 
Or le raisonnement peut s’appliquer a une tranche 
matérielle quelconque. Chaque tranche se trou- 
verait étre ainsi le siége d’une surpression 
moyenne PX?q?/2. On peut done alors consi- 
dérer la pression de radiation comme existant 
tout au long de l’onde sans qu’il y ait lieu d’en- 
visager spécialement un piston moteur vibrant 
ou une surface absorbante quelconque. De cette 
facon, il s’agirait d’une pression au sens ordinaire 
du mot qui s’exercerait aussi bien sur un obstacle 
que sur un élément de la paroi latérale d’un 
tuyau pourvu que cet élément soit lui-méme 
animé de la vitesse u. Cette pression serait alors 
une grandeur scalaire et qui existerait dans tout 
volume élémentaire du milieu matériel au sein 
duquel se fait la propagation et qui se propagerait 
avec la vitesse V. 

Mais, d’autre part, la densité d’énergie de la 
tranche d’abscisse x au repos a pour expression: 


= 
U=+eu2+— Ap [P 
1 


= = 0 X§w* cos*(wt —qx) + 


= PX? q? cos? (w@t—qx), 
U= PX? q? cos*(@t—q~x) 


car P est égal a 0 V*, 0 étant la densité du fluide 
au repos et V étant la vitesse de propagation 
avec w= qV. 

La valeur moyenne de U est = PX3q?. D’ot 
Von conclut a Pégalité de la pression de radiation 
et de la densité d’énergie, d’accord en cela avec 
les résultats admis pour les ondes électromagné- 
tiques. 


b) Mais si les valeurs moyennes de la pression 
de radiation ainsi définie et de la densité de 
Vénergie du faisceau sont égales, il n’en est pas 
de méme de leurs valeurs instantanées puis- 
qu’elles sont régies respectivement par un sin® 
et un cos”. On ne peut vraiment pas dire qu’elles 
soient identiques. 


c) D’ailleurs dans le cas général et en reprenant 
tous les calculs précédents, on obtient successive- 
ment: 
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Par x+xX 
Api=f (: V 
‘Pee x PX) 
=eigrh (: ) V2 ie 
> x BP ” : "” 
d’ou: Ap: = — ya Xf = —off". 


D’autre part, la densité énergétique a pour 
valeur: 


1 1 P2 
ii 2 : PY p ST el 
7 9 of 9P Vp if a of . 


Pour qu'il y ait 4 chaque instant égalité entre les 
deux, il faudrait avoir 


(PY fie 0. 


ce qui ne peut avoir lieu que si f est de la forme 
f=«Vt—x/V, forme a laquelle il est impossible 
de rattacher une fonction sinusoidale. 


d) Cependant on étend le raisonnement aux 
ondes stationnaires et de la facon suivante. Le 
déplacement élémentaire ayant pour expression: 

X= Xp sin(wt+qx)— Xy sin(wt—qx) 


= 2 X)singx coswt, 
la surpression a pour valeur: 
Ap=—P 0X /¢x =—2 PX) qcosqxcoswt. 
L’on a ensuite: 


Ap, =—2 PX,qcosq(x+ X) coswt, 
Ap, =—2 PX,)qcosqxcoswt+ 
+2PX,Xqsingxcoswt+ -, 
Apr=2 PX, X q*?singxcoswt 
=4 PX? q?sin?qx cos? wt. 


De telle sorte que sur un obstacle d’abscisse «= 0 
et sur lequel il y aurait réflexion totale, la pres- 
sion de radiation serait nulle. D’autre part, la 
densité énergétique a pour valeur: 


(Das a -4 X2w* sin? qx sin? wt + 
+ oe 4 PX? q? cos? q x cos? t, 
U=2 PX} q? (sin? qx sin?wt+ cos*q x cos*mt). 


En aucun endroit, U et Ap, ne peuvent étre 
égales 4 tout moment et l’on en est réduit a dire 
que si l’on prend les valeurs moyennes 8 la fois 
vis-a-vis de x et de t, celles-ci sont les mémes. 

Il faudrait ici raisonner séparément sur les 
ondes incidentes et sur les ondes réfléchies mais 
de toutes fagons, il apparait que le raisonnement 
qui a permis d’obtenir Ap, doit pécher par quel- 
que cété. 
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On ne voit d’ailleurs pas pourquoi il faudrait 
toujours considérer des valeurs moyennes alors 
qu’évidemment la pression de radiation a une 
valeur instantanée. 

e) Or, pour établir l’équation de propagation, 
Von a considéré une tranche matérielle limitée 
au repos par deux plans d’abscisses x et x+dx. 
Les déplacements élémentaires étant X et X+dX 
la tranche s’est dilatée de dX et la dilatation 
relative a pour valeur 0=0X/0éx, d’ou la sur- 
pression Ap= POX /éx. 

Mais cette surpression est celle dela tranche dont 
Vabscisse au repos est x et dont l’abscisse actuelle 
est justement x+ X. Il n’y a donc pas lieu d’y 
ajouter Ap,. Chaque tranche matérielle du fluide 
qui oscille de part et d’autre de sa position d’équi- 
libre est donc le siége a chaque instant de la sur- 
pression Ap de valeur moyenne nulle. Et un baro- 
métre minuscule qui oscillerait avec une tranche 
et sans modifier l’état vibratoire du phénoméne 
indiquerait donc une pression moyenne cons- 
tante. 


f) On pourrait peut-étre au contraire se pro- 
poser de calculer la surpression en un point d’abs- 
cisse invariable x, et qui serait mesurable avec 
un barométre fixe. La tranche matérielle du 
fluide qui, 4 un moment donné, est en x, a 
comme abscisse de repos x telle que x+ X= x. 
D’ow pour la surpression correspondante: 


Ap = PX, qcos[wt—q(x,—X)] 
Ap = PX, q[cos(wt—qx,)—qX sin(wt—qx,) +°"]. 
Mais on a d’autre part 
X = Xy)sin[wt—q(x,—X)] 
= X, [sin(@t—q-x,) + qX cos(wt—qx,)+-"], 


eee 
~ 1—qX, cos (wt—qx,) 


= X, sin(wt—qx,) [1+ q X, cos(wt—qx,)+°"]. 
Et en portant dans Ap 
Ap = PXoq[cos(mt—qx,)— 
— q Xp sin?(wt—q-x,)+°]. 


On trouve ainsi qu’en chaque point d’abscisse 
fixe il se produirait non pas une surpression, mais 
une dépression PX? q?sin*(@t—gqx,). La pres- 
sion en un point, du fait méme du mécanisme 
de la propagation serait ainsi la résultante de 
toute une série de termes harmoniques d’ordres 
de plus en plus élevés et qui seraient donnés par 
les termes d’ordres supérieurs, que nous avons 
négligés dans nos développements. 
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2. Mode de calcul proposé par application des 
principes de la dynamique 


Rappelons tout d’abord le mécanisme de la 
propagation d’un ébranlement dans un milieu 
matériel. I] n’est autre que celui que l’on cons- 
tate le long d’un train de chemin de fer. L’énergie 
en se propageant passe alternativement de la 
forme potentielle 4 la forme cinétique. S’il s’agit 
d’un ressort, chaque spire joue a la fois le réle de 
wagon et de ressort tampon. Or ici, aussi petite 
que nous supposions une tranche matérielle, les 
particules élémentaires sont en nombre tel que 
nous ne pouvons pratiquement que raisonner 
simultanément sur un nombre infini de wagons 
élémentaires. Si bien qu’a chaque instant, dans 
cette tranche, l’énergie cinétique est égale a 
lénergie potentielle. Il s’agit bien toujours d’éner- 
gie qui se propage, mais par une course de relais, 
passant d’une particule a la suivante, sans rester 
solidaire de son support matériel comme dans un 
condensateur ordinaire d’énergie cinétique, dont 
Vénergie a pour expression 3 mu*. 

Or une telle énergie, comme toute énergie, est le 
produit de deux facteurs, un facteur de quantité 
et un facteur de qualité. 

Ici ce sont la quantité de mouvement et la 
vitesse. En effet lorsque par exemple un corps 
de vitesse v, arrive sur un corps plus lent de 
vitesse v,, la quantité de mouvement dy se con- 
serve en passant du corps le plus rapide au corps 
le moins rapide et l’énergie mise en jeu dj(v,—v,) 
change de forme. Il en est alors de méme que 
dans la décharge d’un condensateur électrique 
au cours de laquelle énergie mise en jeu a pour 
expression dq (V,—V,) et dont les facteurs sont 
la quantité d’électricité et le potentiel comme 
pour toute énergie électrique. 

Si done, il y a absorption d’ondes progressives 
par un obstacle disposé normalement, l’énergie 
absorbée par unité de surface et en un temps dita 
pour valeur l’énergie comprise dans un cylindre 
de section unitaire et de longueur V dt, c’est-a- 
dire (UV dt), U étant la densité énergétique des 
ondes au contact méme de l’obstacle et V leur 
vitesse de propagation. La variation de vitesse 
de cette énergie est alors V du fait de ’absorption 
et la quantité de mouvement correspondante est 
du=Udt. 

Mais d’autre part, d’aprés la loi fondamentale 
de la dynamique, la force exercée, c’est-a-dire ici 
la pression puisqu’il s’agit de Punité de surface, 
est égale a dyu/dt c’est-a-dire a U, la densité 
d’énergie. Et cela, a chaque instant et quelle que 
soit la nature des ondes incidentes, ondes de 
choe ou ondes sinusoidales et sans que l’on ait a 
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considérer des valeurs moyennes. Et quelles que 
soient aussi les valeurs de la pression et de la sur- 
pression pourvu que l’on prenne pour V la valeur 
de la vitesse de propagation a la hauteur de 
V’obstacle. On sait en effet que cette vitesse croit 
avec elles. 

Dans le cas d’une réflexion totale, la variation 
de la quantité du mouvement étant doublée, la 
pression de radiation est elle-méme doublée. 

Il ressort évidemment de notre facon de rai- 
sonner que la pression de radiation devrait étre 
considérée comme de nature essentiellement vec- 
torielle au méme titre qu’une quantité de mouve- 
ment et qu’elle n’existe et n’est a considérer 
qu’autant qu’il y a obstacle absorbant ou réflé- 
chissant en tout ou en partie. La fagon d’envisa- 
ger la question différe donc totalement de la 
fagon exposée précédemment. 

Si les ondes arrivent obliquement sur l’obstacle, 
il est également trés facile de déterminer l’ex- 
pression de la pression exercée. En effet (Fig. 1), 
la direction de propagation des ondes incidentes 
faisant un angle i avec la normale, la méme quan- 
tité d’énergie au lieu de tomber en un temps dt 
sur une surface s tombe sur une surface S. 


L’action sera donc diminuée dans le rapport s/S 
dans le sens de la propagation et dans le rapport 
s/S+ cost c’est-a-dire cost? dans la direction de 
la normale. Autrement dit, la pression de radia- 
tion sera U cos?i s'il y a absorption, mais en 
méme temps, il y aura une action tangentielle 
égale a (Ucosi-sini). S’il y a réflexion totale, 
les deux actions tangentielles dues respective- 
ment aux ondes incidentes et réfléchies se neutra- 
lisent 4 ’opposé des actions normales qui s’ajou- 
tent. 


3. Calcul par le vecteur radiant 


On sait que dans le cas d’ondes progressives, 
Pénergie qui en un temps dt traverse l’unité de 
surface normale a la direction de propagation a 
pour expression Ap-udt. C’est pourquoi il y a 
lieu de considérer un vecteur dirigé dans le sens 
de la propagation, d’intensité Ap-u et qu’on 
peut appeler vecteur radiant. C’est cette énergie 
Ap:udt que regoit en un temps dt et par unité 
de surface un obstacle placé normalement au fais- 
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ceau. S’il y a absorption la variation de quantité 
de mouvement est Ap-udt/V, d’ou pour la pres- 
sion de radiation l’expression Ap-u/V. 

Il est d’ailleurs facile de voir que dans le cas 
ou nous nous sommes placés, cette expression 
mesure justement la densité d’énergie. On a en 
effet 


Ik , / Ne 9 
Ap:u= yf” et U= 6 [7 Se 


4. Remarques 


Qu’il s’agisse de U ou de Ap-u/V, ces grandeurs 
sont relatives 4 chaque instant a ce qui se passe 
au contact méme de l’obstacle absorbant. Or une 
surface absorbante est mal définie. Une paroi ab- 
sorbante est en effet constituée par des étoffes, 
des tentures, des tapis, des vétements et d’une 
facon générale par des corps alvéolaires. L’absorp- 
tion se fait en profondeur d’une fagon progressive 
et le mécanisme rappelle celui de absorption des 
ondes lumineuses par du noir de fumée. 

Aussi peut-on considérer comme nous l’avons 
fait en somme implicitement que les ondes arri- 
vent sur une surface absorbante animée d’un 
mouvement identique a celui de la tranche maté- 
rielle au sein de laquelle se fait la propagation et 
adjacente a cette surface. Le travail (Ap-u dt) 
exercé par la surpression Ap en un temps dt se 
trouve recueilli par la surface frappée, toute 
Pénergie incidente se trouve ainsi absorbée et il 
n’y a pas réflexion. 

S’il s’agit d’un piston oscillant émetteur d’éner- 
gie, le cas est exactement le méme et le travail du 
piston a justement pour expression (Ap-ud)). 

Dans un cas comme dans l’autre, qu il s’agisse 
d’un phénoméne d’absorption ou d’émission, 
la quantité de mouvement mise en jeu est 
Ap-udt/V et la pression de radiation Ap-u/V. 

S’il s’agit d’ondes stationnaires résultant de la 
réflexion totale d’ondes incidentes arrivant sur 
un obstacle fixe normal a la direction de propa- 
gation, la densité d’énergie au contact de l’obs- 
tacle a pour expression 2 P X?q*?cos*mt. On est 
a un neeud de vibrations et l’énergie se réduit a 
Vénergie de pression résultant de la somme de 
Vénergie de pression des ondes incidentes et de 
celle des ondes réfléchies qui sont égales. Les 
vitesses correspondantes étant de signes con- 
traires et égales, la vitesse résultante est nulle 
ainsi que l’énergie cinétique. 

L’expression Ap-u/V est done nulle puisque u 
Vest. Mais si lon considére séparement les ondes 
incidentes et les ondes réfléchies, on a pour 
chacune d’elles Ap-u/V = P X?q?cos*wt car si 
pour les unes et les autres Ap a la méme valeur 
et si les vitesses u sont égales et de signes con- 
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traires il en est de méme des vitesses V puisque 
les ondes se dirigent en sens inverses. 

On peut donc encore trouver dans le cas de la 
réflexion la valeur de la pression de radiation, a 
Vaide du vecteur radiant; mais A condition de 
considérer séparément les ondes incidentes et les 
ondes réfléchies et d’écrire: 

u 


u = u 


Nous allons voir ci-dessous l’importance de ces 
considérations. 


5. Ondes planes progressant dans un milieu 
absorbant 


On a Vhabitude de ne traiter que le cas relatif 
a des ondes se propageant sans amortissement. 
Or, il n’en est pas toujours ainsi et l’on peut méme 
dire quwil n’en est jamais exactement ainsi. 
Considérons donc le cas d’un milieu absorbant. 
Il s’introduit une force résistante proportionnelle 
en premiére approximation au volume et 4 la 
vitesse; l’équation de propagation est modifiée 
et devient dans le cas le plus simple: 
Ox ox ex 


oe Naame ne teat 


2 


Dans un tuyau il faudrait tenir compte en 
outre de Vinfluence des parois, la modification 
serait plus profonde, il s’ensuivrait par exemple 
une augmentation apparente de @ mais la forme 
générale de l’équation resterait la méme. 

La solution que nous prendrons de forme 
sinusoidale comporte alors un terme exponentiel 
d’amortissement: 


X= X,e-F* + iot—a9), 
Par substitution on a 


—ow*+ajo= P(p+)jq)?= P(?—¢q+ 2jPq) 
2PBhq _2PB 


dot 
o V 


aw=2P6q ou 


a= 


> f2 
soit ae = nz (1 — zl 
q Q q 
a ee i : 
doi y=)" yi—2 V2 “ 
@ q° eas 
= ‘2 ee 
| APS 82?) ° 


Comme on pouvait s’y attendre, la vitesse de 
propagation est diminuée. Mais d’autre part, 
d’autres modifications se produisent. La vitesse u 
des éléments matériels et la surpression ne sont 
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plus en phase, les densités des énergies cinétiques 
et potentielles ne sont plus égales et, en ce qui 
nous intéresse actuellement, l’énergie localisée en 
un point a un moment donné ne se propage pas 
en entier, une partie restant sur place du fait de 
Vabsorption du milieu. 

Et il s’ensuit que U et Ap-u/V ne sont plus ni 
égales ni méme en phase. L’amplitude de U est 
supérieure a celle de expression Ap-u/V. 

Nous ne pouvons donc plus raisonner comme 
au paragraphe [2] pour obtenir la pression de 
radiation: C’est au vecteur radiant que nous 
aurons recours comme au paragraphe [3]. Et nous 
devrons en conclure que la pression de radiation 
exercée par un faisceau arrivant sur un obstacle 
qui l’absorbe est donné par l’expression Ap-u/V. 
Et nous ajouterons que dans le cas particulier ou 
le milieu n’est pas absorbant, cette expression 
mesure la quantité d’énergie au contact de l’obs- 
tacle et que par conséquent cette densité peut 
étre elle-méme prise comme mesure de la pression 
de radiation. Mais il ne s’agit pas la d’un fait 
général. 

Dans le cas owt il y a réflexion, nous devrons 
évidemment tenir compte des considérations du 
paragraphe [4] et considérer séparément la pres- 
sion due aux ondes indicentes et celle due aux 
ondes réfléchies. 


6. Ondes sphériques 


Considérons maintenant le cas des ondes sphé- 
riques. Soit par exemple une sphére pulsante de 
rayon ry trés petit. On pourra écrire pour le poten- 


tiel des vitesses: 
2 


Te 
=A 
3 ie 


qr). 


sin(wt 


D’ou pour la vitesse élémentaire et la surpression: 


oy ra (ees 
uU———_ — A sin(wt—gqr)+ qcos(mt—qr 
Ar 7 |p sin(@t—qr) + qeos(wt—q ) 
=U ,+ Us, 
ain TA 
Ap=o ae 0A—- wecos(wt—qr). 
ot r 


A la surface de la sphére, amplitude de la 
vitesse est A. On a en fait deux ondes qui se con- 
duisent différemment. Pour lune d’elles, u, et Ap 
sont en phase et s’amortissent de la méme fagon 
en fonction de la distance au centre r. 

L’expression Ap-u,/V conserve un signe cons- 
tant et serait égale a la densité d’énergie si cette 
onde était unique. L’amortissement ici n’est pas 
da a une absorption du milieu, mais au fait que le 
volume intéressé pour une méme énergie est d’au- 
tant plus grand que l’on s’éloigne du centre car 
lon a en effet: 
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4ar* Ap-u,=420 A*r} wqcos?(wt—qr) 
= K cos?(wt—qr). 


La valeur moyenne de l’énergie qui passe ainsi : 


a travers une sphére quelconque et par unité de 
temps est constante et égale 4 K/2. Pour l’autre, 
u, et Ap sont en quadrature et si Ap s’amortit 
comme l/r avec la distance, u, le fait comme 1/r’. 
Les conclusions précédentes ne sont plus valables. 
On a comme un certain volume vibrant autour 
de la source puisqu’en un point quelconque 
Ap:u, a une valeur moyenne nulle, l’énergie 
oscillant de part et d’autre du point. On a comme 
un systéme a ondes pseudostationnaires dont 
tous les points ne vibreraient pas en phase. 

Les conclusions relatives a la premiére onde 
considérée ne s’appliquent pas a l’ensemble des 
deux ondes. Le vecteur radiant en un point ne 
sera pas toujours dirigé dans le méme sens et la 
densité d’énergie ne se mesurera pas par l’ex- 
pression Ap-u/V. 

Et c’est cette derniére qu’il faudra encore con- 
sidérer pour avoir la pression de radiation 


4 
=A “3 q cos(wt—qr) X 


Ap:u 
ae 


1 
x E sin (wt—gqr) + qeos(ot—qn)| : 
La valeur moyenne de cette pression est 
1 ke 
= 0A? @ 
Qe pd 


sur un obstacle absorbant a la distance r. Elle 
est égale 4 deux fois la valeur moyenne de la 
densité de l’énergie potentielle. Cette valeur varie 
comme 1|/r?. 

A trés grande distance, on pourra considérer 
Pr comme régie par un cos? en fonction du temps. 

Dans le cas de la réflexion totale, p, est multi- 
pliée par 2. Sinon, la pression due aux ondes 
réfléchies s’obtiendra en multipliant celle due aux 
ondes incidentes par le coefficient de réflexion. 
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Quant a la pression de radiation qui s’exerce 
sur la sphére pulsante elle a pour valeur moyenne 
+ 0 A’r5q*, mais sa valeur instantanée peut étre 
considérable et atteindre la valeur a 0 A*roq. 

C’est dire que la pression exercée sur un obs- 
tacle placé au foyer d’ondes sphériques conver- 


gentes peut étre trés grande. 
(Recu le 7 décembre 1953.) 


Bibliographie 


[1] Brittourn, L., Sur les tensions de radiation. Ann. 
Phys., Paris 4 [1925], 528---586. 
Les tensions de radiation, leur interprétation en mé- 
eanique classique et en relativité. J. Phys. Radium 6 
[1925], 337-353% 
Les pressions et tensions de radiation. Rev. Acoust. 5 
[1936], 99---111. « Les tenseurs en mécanique et en 
élasticité». Masson éd., Paris, 1938. 

[2] Rocarp, Y., La pression de radiation en acoustique. 
Rev. sci. 84 [1946], 329---333. 
«Dynamique générale des vibrations». Masson éd., 
Paris, 1949. 


On pourra également consulter les publications suivantes: 


[3] Baurr, B., Radiation pressure. J. acoust. Soc. Amer. 
19 [1947], 504. 

[4] Borenis, F. E., Acoustic radiation pressure of plane 
compressional waves at oblique incidence. J. acoust. 
Soc. Amer. 24 [1952], 468. 

[5] FuBini-Gurron, E., La tension de radiation acoustique 
et les ondes de grande amplitude. Rev. Acoust. 6 
[1937], 118---127. 

[6] Foxxerr, A. D., Les tensions de radiation acoustique. 
Rev. Acoust. 7 [1938], 73. 

[7] Hertz, G. et Menpe, H., Der Schallstrahlungsdruck 
in Fliissigkeiten. Z. Phys. 114 [1939], 354---367. 

[8] Lucas, R., Sur les pressions de radiation des ondes 
sphériques. C. R. Acad. Sci., Paris 230 [1950], 2004--06. 

[9] MenvoussE, P., On the theory of acoustic radiation 
pressure. Proc. Amer. Acad. Arts Sci. 78 [1950], 
135---184. 

[10] Mercier, J., La pression de radiation en acoustique. 
C. R. Acad. Sci. Paris 232 [1951], 2181---2183; J. Phys. 
Radium 12 [1951], 38S.---40S. 

[11] Ricurer, G., Zur Frage des Schallstrahlungsdruckes. 
Z. Phys. 115 [1940], 97---108. 

[12] ScHarrer, C., Zur Theorie des Schallstrahlungsdruk- 
kes. Ann. Phys. Lpz. 35 [1939], 473---491. 


THE ACCEPTABILITY OF ARTIFICIAL ECHOES 
WITH REVERBERANT SPEECH AND MUSIC 


by A. F. B. Nickson, R.W. Muncey and P. Dusour 


Division of Building Research, Commonwealth Scientific and Industrial Research Organization, Australia. 


Summary 


Groups of subjects seated in a room of short reverberation time (0.15 s) assessed whether or 
not a single artificial echo added to speech or music was disturbing. The artificial echo varied 
50 dB in intensity level and up to 600 ms in delay. The acceptable-echo-level/delay relation is 
shown to consist of three parts; the initial period when the Haas effect is prominent, an interim 
period determined by the decay time of the sound type, and a final period when the level of 
acceptable echo scarcely changes with the delay, being determined by the dynamic range of 
the programme material. 


Sommaire 


Des groupes de sujets assis dans une salle ayant une faible durée de réverbération (0,15 s), 
ont été invités a dire si un écho artificiel simple ajouté a de la parole ou de la musique était ou 
non génant. Cet écho artificiel avait un niveau d’intensité variant jusqu’a 50 dB, et son retard 
atteignait 600 ms. On montre que la relation entre l’intensité d’écho admissible et le retard, se 
compose de trois portions: une portion initiale, ot l’effet Haas est prédominant, une portion inter- 
médiaire déterminée par la durée d’affaiblissement du son employé, et une portion finale ot 
le niveau de l’écho admissible varie 4 peine avec le retard, étant déterminé par le champ du 
dispositif de reproduction. 


Zusammenfassung 


Versuchspersonen, die in einem Raum mit kurzer Nachhallzeit (0,15 s) saBen, stellten fest, 
ob ein einzelnes kiinstliches Echo, das zu Sprache oder Musik hinzugefiigt wurde, stérend wirkte 
oder nicht. Das kiinstliche Echo wurde 50 dB in der Schallintensitét und bis zu 600 ms in der 
Verzogerungszeit geindert. Es wird gezeigt, daf} die Beziehung zwischen der zuliassigen Echo- 
intensitaét und der Verzégerungszeit aus drei Bereichen besteht: dem Anfangsbereich, in welchem 
der Haas-Effekt vorherrscht, einem Zwischenbereich, der durch die Nachhallzeit der Schallart 
bestimmt wird, und einem Endbereich, in welchem die zulissige Echointensitaét sich kaum mit 
der Verzégerungszeit andert, und der durch den Dynamikbereich der Wiedergabe bestimmt wird. 


1. Introduction 


All sounds as commonly heard consist of an 
initial direct sound followed by echoes. Workers 
[1], [2] have studied the acceptability of artifi- 
cially added echoes; the results can be expressed 
graphically as the echo level intensity at various 
delay intervals that are acceptable to a given 
percentage of the audience at the test. Presented 
in this form three parts are apparent: the first 
occupies the initial period to about 50 ms, within 
which period an echo from as loud as the original 
sound to as much as 10 dB greater is acceptable, 
the second occupies the next few hundred ms 
when the tolerability of the echo decreases almost 
linearly with increasing delay, and the third 
when the level of acceptable echo scarcely 
changes with the delay period. 

This report is concerned with listening experi- 
ments using the previous technique of listeners 
seated in a comparatively dead room [2] but the 
sounds used in this study have been recorded in 
reverberant spaces so that in addition to the 


artificial echoes there are natural echoes in the 
sound. The results lead to a hypothesis which 
explains them in terms of the characteristics of 
the original sound. 


2. Experimental technique 


The equipment used in this work is essentially 
that used previously [2]. The endless tape system 
is unchanged but the amplifiers and the control 
system have been reconstructed to provide in- 
creased operating facility and somewhat improv- 
ed sound crispness and quality. A new room for 
listening has been constructed; it measures about 
6m <x 5m x 3m and has a reverberation time/ 
frequency characteristic as shown in Fig. 1, 
room A. Two loudspeakers placed side by side 
were used for the initial and echo sounds and the 
number of listeners who answered the question 
‘*‘are you disturbed by the echo ?” in the affirma- 
tive was recorded electrically as before. 

Several sound types were used; three were one 
excerpt of 20s length read at 5 syllables per 
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Fig. 1. Room reverberation times. 


second in the three rooms numbered A, B and 
C in Fig. 1, the fourth was string music recorded 
from a quartette playing outdoors, the fifth this 
string music played in room C and recorded from 
amicrophone, and a sixth was the speech recorded 
in room A re-recorded concurrently with random 
noise onto the endless tape. As results show, it 
is unlikely that the actual words and music are 
significant; they are therefore not reported. 

The audience at each listening test was nor- 
mally ten; occasionally the whole ten were not 
available and the corresponding adjustment was 
made in the evaluation of the results. The eight 
persons of the acoustics group of this Division 
were included in the ten and, from the previous 
paper, it can be expected that their results will 
be representative generally. 

The sound level used was the same as in the 
earlier test, being about 80 dB (peak) above 
0.0002 dynes/cm?, as desired by the subjects for 
comfortable listening. In addition, for the speech 
from room A a further test was conducted at a 
level 20 dB less to find the effect of intensity on 
the results. 

Finally, the overall level of the non-echoed 
sound of each type was recorded by a Briel and 
Kjer type 2301 high-speed level recorder ad- 
justed to have an instrument decay equivalent 
to a reverberation time of 0.14 s, being close to 
that of the ear [3]. The rate of sound decay 
after each syllable of speech or music has been 
termed the programme reverberation time. 


3. Results 


The results obtained from the various sound 
types are shown in Figs. 2 and 3. Initially a plot 
was made of the percentage disturbed by the 
echo against the amount of the delay and from 
these plots the echo level/delay relation accept- 
able to 20 and 50 percent of listeners was deter- 
mined. The precision of the results decreased 
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with increasing reverberation time on the original 
sound and was greater for speech than for music. 
The standard error for a single observation for 
speech recorded in room B was about 12 per cent 
for any particular level/delay combination. Meas- 
urements were repeated to reduce the error of 
the mean to about 5 per cent for speech and 10 
per cent for music. 

The records of the sound level were analysed 
to measure the level difference between the peaks 
and the troughs (roughly within one reverbera- 
tion time) and the minimum noise level and also 
the programme reverberation time measured as 
each syllable or sentence of speech or music died 
away. These are recorded in Table I together 
with details of the equivalent reverberation time 
(obtained as a tangent to the echo level/delay 
time curve), and the level below which echo 
acceptability is independent of delay taken from 
Figs. 2 and 3. 


10 ;= ] ail 
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Fig. 2. Acceptable echo levels for speech. 
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Fig. 3. Acceptable echo levels for string music. 
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As might be excepted, the result of the test at 
a level 20 dB below that normally used gave 
results which did not differ essentially from those 
at the usual level. 


Table I 


Listening test 
results 


(20% disturbed) 


Level records 


Equi- Level for 

Sound type valent echo 
reverbe-| peaks reverbe- | indepen- 

ration and ration dent 
time | troughs time of delay 

s dB dB 


Speech 

room A 0.20 42 42 0.18 —45 
Speech 

room A with 

added noise 0.20 28 28 0.18 —35 
Speech 

room B 0.5 35 35 0.52 —30 
Speech 

room C dt 32 32 1.05 —30 
String music 

outdoors 0.35 23 38 0.3 —25 
String music 

room C ibsul 20 32 0.95 —20 


4. Discussion 


It is quite apparent that the nature of the 
original sound has a profound effect on the 
results obtained in listening tests such as de- 
scribed here. The agreement of the programme 
reverberation times and maximum intensity chan- 
ges, measured from level records of the sound, 
with the equivalent reverberation times and mini- 
mum echoes corresponding to the more critical 
section of listeners (20 percent disturbed) is quite 
striking. In fact it appears almost certain that 
it is the former two parameters that determine 
the shape of the curve after the initial period 
when the “Haas” effect operates i. e. when the 
ear is capable of accepting echoes greater than 
those given by a linear reverberation decay. The 
tentative criterion for auditoria proposed by 
Bor and Doak [4] follows this idea. 

As a further check on this hypothesis the 
recorded sounds used in the previous study [2] 
have been examined and the results are compared 
in Table II. Again the agreement is very good 
with the exception of the slow music. For all 
cases the ratio between peaks and troughs of the 
sound was measured also after passing through 
a tuned circuit: for the slow string music the 
ratio was then about 28dB and for the slow 
organ music 20 dB, whereas it was practically 


A. F. B. NICKSON, R. W. MUNCEY and P. DUBOUT: ACCEPTABILITY OF ECHOES 


449 


unchanged for the other sounds. This suggests 
that the frequency sensitivity of the ear is being 
used to detect low level echoes in slow music. 
Further, many of the subjects complained of the 
difficulty of detecting echoes with slow music 
and about 30 percent reported being disturbed 
by an echo even when the sound was presented 
without echo. These results would appear to be 
less reliable than the others. 

An independent confirmation is supplied by 
MEYER and ScHODDER [5]. Study of their Fig. 8 
shows that for echoes to pass unnoticed by the 
more critical section of an audience they must 
fall close to the level expected from a linear decay 
at a rate corresponding to the reverberation time 
of the listening room. 


Table II 


Level records 


Listening test 


Pro- Ratio Equi- Level for 
Sound type gramme between valent echo 
reverbera- | peaks and | reverbera- | independent 
tion time troughs tion time of delay 


Speech 


fast 0.20 —35 

slow 0.28 —34 
String music 

fast 0.35 —32 

slow 0.4 —30 
Organ music 

fast 0.6 —31 

slow 0.8 —30 


The results from listening in a dead space to a 
sound recorded in a live space are qualitatively 
comparable with those previously reported for 
listening in a live space to material recorded in 
a dead space. For example, replotting Haas’ 
Fig. 4 in the form used here suggests (from the 
slope of the tangent) an equivalent reverberation 
time of about 0.8 s. In both instances then, the 
tolerance for echoes is largely defined by the 
masking echoes introduced by the room. This 
leads to the conclusion that, apart from short 
delay echoes, the listening tests are only a com- 
plex method of measuring phenomena connected 
with reverberation. 


5. Conclusions 


The relation between acceptable echo level and 
echo delay in artificial echo listening tests de- 
pends in the middle period on the natural echoes 
present in the recording room, and in the long 
delay period on the difference in sound level 
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between the loud and quiet portions of the pro- 
gramme material. Any attempt to translate 
results of listening tests into criteria for examina- 


tion of pulse measurements from rooms would’ 


therefore produce merely a reflection of the test 
conditions and not be universally accurate. 
(Received 15th December, 1953.) 
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LETTER TO THE EDITOR 


concerning the article by 
“M.-L. EXNER and W.HAmPE, Experimental determination of the damping of pulsating air bubbles in water” 
(Acustica 3 [1953], 67) 


by M. STRASBERG 
David Taylor Model Basin, Washington 7, D.C., USA. 


M.-L. ExNER and W. HAmps, in the above mentioned 
article, report the existence of ‘“‘anomalous” bubbles. This 
letter discusses a possible explanation of the phenomenon. 

An anomalous bubble is characterized by the following 
acoustic behaviour: 

1. The measured resonant frequency does not agree 
with the frequency calculated from the measured diameter 
of the bubble using MINNAERT’S relation between frequency 
and diameter. 

2. The damping at the measured resonant frequency is 
abnormally low. 

3. The measured resonant frequency may suddenly 
change. 

4, The bubble usually has a dust particle on its surface. 

The article referred to examines as possible explanations 
of this behaviour the influence of the dust particles on the 
mass reactance of the water, and on the surface tension of 
the bubble surface, but concludes that neither of these 
explanations is valid. As an alternate tentative explana- 
tion, I suggest that the behaviour may be associated with 
surface oscillations of the bubble. 

Surface oscillations of a bubble are usually of such low 
frequency that they are neglected in acoustics. However, 
for very small bubbles, the frequency of surface oscillation 
may be of the same order as the frequency of ordinary 
volume pulsation. The circular frequency ~,, of the surface 
oscillation of a bubble of gas in a liquid is given by LAms [3] 


. 
, w2 = (n2—1) (n+ 2) T/oR’, 


where n is the order of the surface harmonic describing the 
shape of the oscillating surface (n =2 is the lowest order 
possible for a surface oscillation), T is the surface tension, 
is the liquid density and R the mean radius of the bubble. 
For comparison, the frequency , of volume pulsation of 
a bubble is given by MINNAERT [2] as 


wo} = 3y Plo RY, 


where P is the mean pressure and y the ratio of specific 
heats of the gas inside the bubble. The ratio of the surface- 
oscillation frequency to volume-pulsation frequency is thus 


(®n/o)? = (n?—1) (n+ 2) T/3 y PR. 


It is apparent that as the radius R becomes small the sur- 
face-oscillation frequency approaches, and may indeed 
exceed, the volume-pulsation frequency. 

As an example, consider one specific bubble described 
by Exner and Hamper which, at a radius of 5.5x 10 em, 
had a resonant frequency which suddenly changed from 
about 81ke/s to 97ke/s. These frequencies correspond 
approximately to the ninth- and tenth-order surface modes. 
The sudden change in observed frequency could be a result 
of a change in excitation from one mode to another. 

The reduced damping of these anomalous bubbles may 
be due to the smaller radiation damping associated with 
surface oscillation. The damping due to heat conduction 
would also be smaller; in fact, probably the only important 
source of damping is viscosity. 

The excitation of these surface oscillations by a uniform 
sound field may require some non-symmetry. Non-symme- 
try may be supplied by the dust particles attached to the 
anomalous bubbles. ; 

A test of this explanation of the anomalous bubbles 
could be performed using the apparatus of EXNER and 
Hampg, by determining whether a single anomalous bubble 
could be excited at a series of frequencies distributed 
according to LAmB’s formula given above. 

(Received 2nd May, 1954.) 


[The English editor comments that he has observed 
flutings on the sides of the large cavities formed at impact 
of a sphere on a water surface which he ascribes to higher 
modes of these vibrations given by Lamp’s formula. 
KorNFELD and Suvorov (J. appl. Phys. 5 [1944], 495) 
observed similar effects on small cavitation bubbles. Pre- 
sumably these higher surface modes involving no change 
in volume are inefficient radiators of sound.] 
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